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本 书 向 读者 介绍 三 种 数学 工具 ， 即 变 分 法 、 有 限 元 法 和 外 推 
法 。 它 们 在 自然 科学 和 工程 技术 中 有 广泛 的 应 用 。 

变 分 法 具有 理论 的 完整 性 和 应 用 的 普遍 福 。 它 用 统一 的 变 分 
原理 解释 了 无 数 的 自然 现象 ， 又 用 巧妙 的 计算 解决 了 多 种 实际 问 
题 。 无 论 是 在 理论 力学 、 弹 性 力学 和 塑性 力学 中 ， 还 是 在 电磁 场 
理论 、 相 对 论 和 量子 力学 中 ， 只 要 应 用 变 分 原理 ， 就 使 这 些 学 科 
显示 上 出 内 在 的 和 谐 。 在 本 书 第 二 章 中 ， 系 统 地 秀 述 了 古典 变 分 学 
的 基本 内 容 。 读 完了 这 一 章 ， 您 就 可 以 用 变 分 法 来 解决 工程 技术 
问题 了 。 

有 限 元 法 是 最 近 二 十 多 年 来 最 引 人 注 意 的 计算 方法 。 有 趣 的 
是 ， 这 个 重要 方法 是 分 别 由 数学 家 、 物 理学 家 和 工程 师 独 立地 创 
立 的 。 也 就 是 说 ， 它 有 三 个 来 源 。 目 前 在 我 国务 大 学 的 非 数 学 专 
业 中 讲 的 有 限 元 潜 ， 几 乎 都 是 从 沫 、 板 等 具体 问题 引入 算法 ， 着 
重 于 实例 。 这 样 使 学 生 能 较 快 地 学 会 有 限 元 算法 。 但 是 ， 这 样 佑 
也 有 一 些 缺 陷 ， 使 毕业 生 很 难 在 有 上限 元 理论 方面 进行 研究 ， 对 有 
腿 元 结果 的 误差 合计 也 很 难 掌握 。 对 当代 正在 发 展 的 各 种 有 限 元 
方法 如 非 协 调 元 、 杂 交 元 、 边 界 元 ……》 较 难 接 受 。 为 弥补 这 
一 缺陷 ， 在 本 书 第 四 章 中 主要 是 讲 有 限 元 法 的 数学 思想 ， 着 重 指 
出 有 限 元 法 是 在 变 分 法 这 棵 二 树枝 头 开 放 的 新 花 ! 看 了 这 一 章 ， 
对 您 今后 攻读 有 限 元 理论 也 许 是 有 好 处 的 。 

外 推 法 是 一 种 很 好 的 近似 计算 方法 。 对 于 已 求 得 的 低 精 度 近 
似 信 ， 只 要 作 儿 次 最 简单 的 四 则 运算 ， 立 即 得 到 高 精度 的 近似 
È 可 以 说 它 是 对 近 级 值 进行 “ 精 加 工 ” 的 妙法 ! 就 是 在 广泛 使 
用 电子 计算 机 的 国家 中 ， 技 术 人 员 对 外 推 法 也 十 分 欢迎 。 因 为 这 
种 算法 可 以 大 大 缩短 计算 时间， 程序 X $5. PE 1910 年 ， 
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L.F.Richardson 就 提出 了 外 推 法 ， 但 多 年 来 它 几 乎 只 有 时 来 解 
次 一 维 问 题 。1965 年 苏联 学 者 在 侦 徽 分 方程 的 外 推 方面 有 了 
ZEEE. 19834, 3A X A B AROE XE WD UH Eden] (XEEREK 
外 推 》 一 书 的 英文 本 在 西方 出 版 ， 此 书 在 局 美 数学 界 引 起 震动 。 
外 排 法 的 研究 又 活跃 起 来 。1980 秆 以 来 ， 中 国学 者 林 和 群 、 陈 竺 
葬 、 吕 洲 、 沈 树 民 等 人 研究 有 了 腿 元 解 的 外 推 ， 取 得 巨大 进展 ， 写 
下 了 外 推 法 中 的 新 篇 章 ， 引 起 国内 外 的 重视 。 但 是 在 国内 的 教科 
书 中 ， 半 于 偏 微分 方程 数值 解 的 外 推 法 还 是 空白 区 。 

在 本 书 第 五 童 中 介绍 了 外 推 法 种 它 在 仿 微 分 方程 中 的 应 用， 
试图 填补 这 一 空白 。 

只 要 具有 工学 院 本 科 二 年 级 的 数学 程度 ， 就 可 以 看 懂 这 本 
书 。 如 果 您 只 想 学 外 推 活 ， 可 直接 看 第 五 章 。 如 果 您 关心 的 是 有 
限 元 法 的 数学 思想 ， 可 只 看 第 一 、 四 章 。 要 尾 您 对 变 分 法 的 应 用 
感 兴 趣 ， 可 只 看 第 一 、 二 ， 三 章 。 不 过 ， 我 还 是 劲 您 在 看 第 五 章 
外 推 法 ， 它 对 您 的 业务 工作 也 许 很 有 用 。 


作者 
19844F12 H, 
于 华东 交通 大 学 。 
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第 一 部 分 “ 变 分 法 和 有 限 元 法 


第 一 章 5l iE 


Aci — ECIEEA IE Ha un 

ERU E EAS- T EOE RERBA. PUE T BD 
VNB , KBER: KA L EE Ak E e OM GR £u 5E 
转 ， 问 桶 中 水 面 呈 入 么 形状 时 所 受阻 力 最 小 ? 他 把 这 个 向 题 归结 
为 : 求 函 数 y 一 帮工 ) 


Whi MECA LANDE A, 但 牛顿 当年 用 来 解决 
河 题 的 方法 还 不 是 变 分 学 中 的 基本 方法 。 

导致 变 分 法 建立 的 著名 问题 是 瑞士 数学 家 约翰 、 贝 努 里 在 
1696 年 提出 的 ， 即 “最 速 降 线 问 题 ”。 此 问题 发 表 于 当年 6 HE 
的 和 教师 学 报 少 上 。 问 题 一 提出 来 ， 立 即 吸 引 了 当时 世界 上 许多 
最 卓越 的 数学 家 的 注意 。 在 该 杂志 1697 年 5 月 号 上 刊登 了 牛顿 、 
莱 布 尼 慈 、 党 比 达 和 贝 努 里 兄弟 名 的 解法 。 他 们 殊途同归 ， 轴 各 
不 相同 的 方法 得 到 问 一 解答。 管 案 是 有 趣 的 ， 原 来 “最 速 降 线 ” 
居然 是 摆 线 ! 这 些 解 题 方 法 中 列 含 着 天 才 的 思想 ， 就 像 负 矿 中 储 
芒 肴 铀 -235 一 样 。 但 是 只 靠 这些 方 法 ， 变 分 法 是 创立 不 起 来 的 。 
正如 只 依靠 古代 的 极限 思想 无 法 建立 微分 法 一 祥 。 那 么 ， 谁 来 滞 
除 杂质 ， 炼 出 真 金 呢 ? 历史 安排 了 天 数学 家 尤 拉 1 他 是 创造 各 种 
数学 方法 的 大师， 他 创立 的 方法 几乎 总 带 有 明晰 性 和 一 般 社 。 他 
在 1784 年 解决 了 更 广泛 的 最 速 降 线 问 题 ， 但 他 对 自己 当时 用 的 方 
法 还 不 满意 。 他 开始 寻找 解决 这 类 问题 的 一 般 方 法 ， 这 种 方法 终 


D *JUsmÓu€"iBlme Bernoulli 和 好 前 弟弟 John Bernoulli, 
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于 被 从 找到 了 ， 蛮 分 潜 也 就 建立 了 。 无 拉 在 1736 年 的 论文 中 指 

Hi. BEH 


Qcy)=) FC Go, y! G2342 (1) 
取 家 人 入， 函数 y+) 必须 满足 
F df 9F Y. » 
A O O c i 


方程 《2) WERA ERNI” o P, MEERKATTE 
的 第 一 个 证 明 竹 往 失 之 繁琐 ， 睁 拉 这 次 也 不 鲍 外 。 后 来 法 国 太 数 
学 家 拉 格 兰 日 改进 了 尤 拉 的 证 明 ， 人 性 之 完全 符合 数学 分 析 的 精神 
而 且 二 分 简洁 ， 他 在 1755 年 把 这 个 证 明 告 激 了 尤 拉 。 至 此 ， 变 分 
学 作为 一 个 新 的 数学 分 支 算是 形成 了 。 

刚才 谈 到 的 是 变 分 法 古典 部 分 的 内 容 ， 它 所 研 避 的 主要 问题 
可 以 归结 为 ， 在 适当 的 函数 集合 内 选择 函数 y=/(7*)， 使 积分 
C1》 取 极 值 。 解 决 这 一 问题 义 归 结 为 解 苑 拉 方程 2》。 看 起 
来 河 题 似乎 并 不 复杂 ， 办 法 也 很 平常 ， 不 过 和解 一 个 微分 方程 罢 
了 。 其实 不 然 ， 我 们 依靠 这 种 方法 ， 就 能 用 统一 的 数学 程序 来 解 
决 自 然 界 和 其 他 方面 的 千 获 万 吊 的 问题 ， 就 能 用 奇妙 的 变 分 原型 
解释 无 数 的 自然 现象。 


£ 考 x E 
(1: M- XB, Chr Er dufR. PAFA 1980, 
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第 二 章 ” 变 分 法 初步 


在 这 里 不 介绍 泛 函 的 精确 定义 ， 因 为 那 祥 做 震 站 比较 多 的 现 
代数 学 知识 ， 势 必 兜 大 轿子 。 实 际 上 ， 用 几 个 例子 出 可 以 说 明 
“ 泛 亲 ”这 一 概念 的 实质 。 

例 1 HARY- EAEL Oo ,11 上 连续 ， 且 f(z 
0 ， 风 微分 


QUf3- [ /ra 


是 曲线 y= f(xi5iEESu-OI[RBIBOEID, GR, o FE A 
fELO,1JLiEBEHSESBPMEUu-fGoXi, QC TB RE Mia 
mie. Ir c5 


fix)mm 
jul 
. n t 
Wi — eip- 
Ot fir) NES 
ES fix)—-V r 
1 二 
Wl QC Gy» | v5 dz =2 


MERRET TERE HATERA CAU. 
r ERE Du (Ru. VERA RIRI, WE ERDER 
值 ， 2 部 有 淮 一 确定 的 值 与 它 对 上 应， 那么 我 们 就 称 变 量 = 为 口 变 
H. Anivy WATE LABELS Ej DOE XE UB H H EH x 


Q o Ffig fih 5s, O f tsp ky- fial, Hayn, Xi 
r-aE 'j5ir-5gsztmE H EEA 
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e doc: 
EE WEE. RZE OACACSCNÉBU. HOEWilt CUuE— Wü 
FATE. PEWA 1, SERTO ET, QLI Cr RE UE 
HE” T, KAREHA? 这 说 朋 我 们 可 以 把 变量 QL] 
看 成 是 “ 丽 数 六 了) 的 函数 ”。 这 人 [六 2 也 记 为 GCC 门 ， 在 此 ， 
f3xunzeB, QE f cu. 
$012 我们 把 一 切 在 [ea，B] 上 连续 的 图 数 GE X BE XE. LÀ 
TiO HEERA, bL FERS 


MT J= | feoda 


ER ^ TELa,52. EXE SER p d o MORTE A EAS E 
HRS MURA FRZ EEA 1f. f BAIERI S5. 
我 说 丰 出 在 邓 和 上 了 之 间 ， 有 一 种 摧 数 关系 ， 提 这 是 推广 了 的 
函数 关系 ， 也 就 是 下 文中 要 讲 的 “ 泛 丽 2 。 诛 来 所 谓 “ 兴 盟 ”， 
不 过 是 更 广泛 意 交 下 的 画 数 关 系 巡 了 。 
例 5 把 fx) 在 Ca， bI ERSTER EDIBIENCADS, MER 
一 个 在 Ca，6561 上 有 定义 的 实 落 数 /(z) ,六 有 确定 的 值 与 它 对 应 。 
AER, JACO, 2x2, W 
NÜcosr^— 2 
NC Csinr]— 3 
N(x*-11— 


从 以 上 风 个 例子 中 ， 我 们 可 以 概括 出 如 下 的 重要 概念 : 把 具 
ie MR PERS AU RO TED, TEA D:pmBHERB EIU 
《 即 对 任何 f ED》 ， 变 量 久 都 有 了 唯一 确定 的 值 与 它 对 应， 肝 么 
FE Q ui ir FoU PEN CEMTHZ US. UA 

Q-—QUf(x)) 或 Q=QL/JƏ 
3X dé iz BREED T TUE RES. HATRA RKK, IX PH E 
经 够 用 了 。 


Q Pimi pham tE, mii = a-p, DENTES. 
© QrofconitiifgEQUsuix d, 在 此 WW = y(r) = fiz), 
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1. iXi&QryCGOn m. uCr)oB HE 

对 于 泛 画 日 [yx)1,ytX) 是 集合 司 中 任何 元 素 。 召 果 y(3) 册 
yox BE uy CE) Myr) gy COR] Ry e EYT) 上 的 变 分 
(Gy (02 y Cx HS 3 DO ， 记 作 


àyg-—yiz)—1 o) C32 
我 们 以 后 常用 u 
óycy(cx)—y) C42 
这 是 指 在 y(+) 上 的 变 分 。 
2. HRPA 


X PiEHOLyCOONE, WEZYR EASy 充分 小 时 ， 
0 的 改变 量 可 以 任意 小 ， 那 么 就 称 泛 画 QLy(z)] 是 连续 的 。 例 
Qty G3 [ God 


Tir. Hyl) ECL, b) (SAFE S 151250 Qly] 
有 定 尺 。 对 性 给 的 £2>0， 只 各 


sry (x i 
Ed yo [LÀ 


Cury RT cLa, OD 


n 
| QC G3 7 QCy C2 i EG yr) dz 


«f toco -ucnldz «tans 


可 见 QceGone | uo daR EREE, 


3. fETEIT PH 
ipii au Q[y(%)J 与 y(z) 的 关系 是 线性 的 ， 也 就 是 Wi. OQ 


蕉 动 论坛 __ _ ~ ~ ~ htp 人 vibbiceducnivipbbs 
+ 6 " 
识 足 以 下 条 件 
120 Qrey(x))-2CO[uyGr)] ec 为 手 交 常数? 
2) Qty, GO FE gr) — QUy Ce) QU (Gr) D 
XB, XkQryCOr)D*W£ETEI;EPS, [dn 
Qtyc)n- yeda 
ERTZA, PRE 


Olele) terlay m Cey Ge ey GOMz 


me wGMdzt eid unte e (QC GH) 


te; QUyG 3 
但 是 


Qty [God 


不 是 线 件 泛 喇 ， 这 一 点 请 读者 自 己 验 证 。 

4. TEE 

EEEH "r9" ERASO FBUESUSTEIO, MER 
PEZH., ME ANED” EE RAPA” XX— 8pm 
的 推广 。 如 果 明 和 白 了 这 一 点 ， 您 名 懂 得 了 EI Wapi. 
为 了 把 间 题 说 和 清楚， 我 们 有 必要 先 复 习 一 下 “ 短 分 ”这 一 报 念 。 

现在 我 扣 一 个 看 来 非常 简单 的 问题 (PA Ec CO 的 
微分 ? 

IREA tpTS. kA LAE 

du-—f'irr 

虽然 我 不 能 说 您 答 诺 了 ， 但 是 您 并 没有 说 清楚 “微分 ”的 实质 。 
要 是 您 在 谈 到 微分 定 闷 时， 只 能 说 出 上 天 这 徊 话 来 ， 所 各 嘲 我 自 
言 ， 您 并 没有 学 好 数学 分 析 ， 

究竟 什么 是 函数 的 微分 昵 ? 我 们 知道 阴 数 关系 y= 二 了 (xz) 往往 


Q uxlsHPum—GXTGEeGzik., HI 
Ecpat) +t cay. 38 eu Cy Ux) Dt cu QU GO ZR EE (DER 01,0. 成 立 。 


36. zh 16-18. http://vib.hit.edu.cn/vibbbs 
e 7 a 
Rd EAS. HERAA COORMEAy=fle+ Ar) xr) 
EEE taaa. mi. ~ 


上 


EMIT SE oO Bs 
E a AX ， ÅT 
| | : ， 


Å YSN ÅT) -Snr 2 &in 


"rc Miba ^ reai Aut A Ayi Ar qpE 
or -JERR ERU. EAE Ar a R: DR Ms =f) 
定点 工 处 有 导数 广 (>)， 这 时 


din we EE 
TE 
AE —f':x)kbe  QUuXr- 0t, a 0) 
Hr 以 


Ay—jf'(x)Ax-—-aAr-dy-Jda.Ax (8) 
TARERE, STejpAGUES ERR B PE Bm 
dy- dis ar MEGER CORR SOG Mudy Ae nk 
例 芍 。 通 常 说 dy EA p eR EL, 352 34ba Ax, 4 år- 时 
4-0. Br 

-VAT 0 ChAr 0) 
"p oAzxdkpEA.HERÉI ZUM. UB 

EA går} (ATO 

Wr Er 

Ay= F v EAS PT) (6) 
WD dyc-f'GOo^c e Eo dtu E Spr LER DLE (a 
既是 Ax# 的 线性 部 分 KEAVAT EA MAA “REEE ND 
P" o. XT UE REOR CA e 
Sr. IAFB A BEBIT dep) EC. 

HEH A. "Ax AADA E RGAN JB HA dy 
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T (7z)Az 作 为 增 量 Az 的 近似 值 了 。 再 来 看 = Simnz 这 个 函数 ， 就 
有 起 下 近似 式 
Ag —sin(r4 AÀx)-—sinref'(rjJAx-coSr-:AT 

用 cos zx。AX 作 为 Ay 的 近似 值 ， 这 就 方便 多 了 。 可 以 说 ， 用 线性 
关系 来 逼近 非 线 性 关系 ， 这 是 微 积分 学 的 基本 思想 1 

有 了 以 上 的 锥 备 ， 我 们 就 可 以 讲 认 图 的 变 分 。 先 看 一 例 ， 
APEM QFy(z)] 一人 oz(z)dz， 它 的 增 量 可 表 为 

AQ-—Qty,(z)1—- QUyCGxO)1 Lyle) t öy] QUuCx)1 


=f Cytr) +y dz -È yoda 
=È Lya) +2yle dy * (öy de f yde 


=f 2gC(x) óydz e ayaz 


年 此 åy =y r) yir) 
uu, EADARRA MSM. A 


| 2yr jðydz =T) y] 
pd Ey BER. Tiylr) dy l EA Tóva S BA. a E 
T(y (G2), Côyl=| 2yCz)Casdz 
=C ['2yóydz CTty (63 
ii 
Tuir, 8y t óy. =È 2y EXS y toy dz 
zi ayCz) oy dz ['296089idz 
=F {ylz} ðu] + TlyCG,6g.J 
们 再 来 观察 第 二 项 | .(6y)2dz。 在 此 63 一 aiKz) 一 5(z)， 其 中 
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yCx)E DABEI. vi OO RECURIS ENS vO0OE i) EIS 
TCLa,53 
若 max ly,Gr)—gy(r) —maxlóy| =-— 0 


gc 


zx max (ðy) (b-a) 
b 


QUERI 


| Cayytds 


"[4n 
ME 


~ max lèy] | 


——— D 


这 就 是 说 ， 当 max fyi -— 0 Hf. [5w)*dz 是 出 max Sy! 高 价 
HEADE, ARES RUA 
| G8).e= OLAN) 
TE 
AQ =T Cr) dyl Sy) 
这 个 公式 与 《5 ) 式 何 等 租 似 1 e a Ae a OX, BU i2 PA 


Q tyco1- [ f Godz 的 增 量 可 分 解 为 两 部 分 ， 第 一 部 分 是 


6 的 线性 泛 芍 ， 第 二 部 分 是 比 60 更 高 阶 的 无 穷 小 量 。 这 种 情 形 
ipe AGE RES AA GL GEI 征 情 的 读者 立 冯 会 想到 ，“ 这 第 一 
MATEY) AVEREA HEHn? 正 是 这 样 1 现在 我 们 
WERDE “ARH TE” EEN, 

EX XyiBO[(syOs) u OARE Sy CHIC OR) 
分 ， 参 看 本 章 身 2 中 的 第 C2 )》 SO. "izIQO RGB 5QO-Q 
[yCz)d óy OLDO h ABL AORn pg 

AQ = Tyer dyl BUg(x),697 (7) 

Ee, TÜy(x).0g2XJórdn ChvCOHB$EO 是 线性 泛 函 而 


BLyCx) 0v] , 0 Chy(Cc)[BhigE, max lóyl ——0) 
max itl 


那么 ，7 [TY AOR MZE E WIES., WIE Qv 
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(Ope sQEORM med €ipkdgDEGM SU. m CEGEBOR 
CIERRA E RREI MA AA E A ERA E A 
PET A EE ERR HEU DIE, 
再 举 一 例 。 将 在 [0 . 1] tA ELTE ERR FE x 记 ED. 1L 
Myn ED, yt 0 ) 的 信和 是 唯一 确定 的 ， 饭 以 Qty(2)]=y? €0) 
是 定义 在 口上 的 一 个 汉 函 。 当 wt7) 取 定 ， 求 此 泛 孙 但 的 变 分 60Q。 
显然 
AQz-Qry(r)--6y]— Onr -[y(CzY- áp 
— y^(0)- yO - 2y(0Y-8pg(C0) 1 ógyCO).T 
— y*(05 -2pt00óyuCO Ff 8g(002* 
TI Sy = dy hees SEET — YET) 2... 
3238  AQ-2y(008g0) Ay W AR HI 
29(0)8gC0) — TÉ yr 6p HEX óy B 2X YE iE E. LÀ 
max jðy] —— 0 时 ， 
(Sy(0))  Cmax fóy Y 
max|óy] ^ max [ow 
FEMA MERAQ AREER REOR dE. BD 
GQ =2y 640) 
E, WTA =y, B6 
du-2y(r)dr 
此 微分 式 和 上 面 的 变 分 式 是 多 么 相似 RERA SES RSS 
与 栈 数 的 导数 也 许 在 计算 方面 有 联系 。 惠 展 果 然 如 此 。 现 在 来 介 
mR, CHEE SRO) 的 变 分 存在 ， 
xem eidyEEvCOmóyg-u0o-vy M WR 
Qrgy(r)raón5) 
显然 ， 它 是 变量 o 的 函数 ， 可 记 为 
QUu(Gr)c-aóg'-vc0D 
REME, PI HEP eCGOTEe- 0 处 的 导数 和 00) 就 等 于 泛 国 
QtuCr)JEg2E 4, TüYESEUEH IE BR. 
=E WAEL AE l Tyr) Sy E 在 ， 


-max lAn] —20 
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BE 24 UE SEARAET OM ERO = Oruro -edy dA ERE 04b 
BH (yCGOBOUSNBED. 
W3H. IaisgCo) ENS r QUsG)TGÓvy]-v(a) 
wap K, H 


ataim an 
sry G 
-qim QO r eriy O 
rz —3 il a 


HAQ yG NA Ee, Mre 
QLyCGr)-óygi—QLyCr] 
Tytr), on. + Pyl) öy] 
于 是 
QUyGe)- ady! GL yr)) 
-T(yGOo.aày)4 Bi gCx),aódy) 
Hy ME, TüyOr)àsgixXósminiihfEited. Bre 
TIytr) vou =af Ly(r), Oy) 
lim QCa)— eto) 


a —M 0 G 


— jm ZUe dt BUyCr),aóy] 


OO 


z- f a 


T yiri yit D ylz ay 
DE yir) SEE P YRT) GOY 
Do Ü 


—fiy(r),8gy)e lim 
T 一 名 [D a 


但 是 ， 册 公式 《7)》 可 知 


Pey yl 0 Quy, max [Bub 一 >0) 
max ló yl 


Plyler) aòn] BÜgCr),aóy, | max lad 


lim -lim [AE 

a- ü F " cU max iul a 
E max lsu 三 T t lim B:yix UAM. 
d [-mas nau ; 7 max lósi Of, Hg | 0 
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sp LO J=T ytz Ay] (89) 
5. RAH 
BiEBQUyCGOMEBBZuGOm fc EREM, Biyo 均 
位 于 平面 上 上 的 有 界 区 域 号 内 【不 妨 认 
为 刀 是 矩形 ) ， 如 图 1 Bp. 


MWER QOL) EBD S 
Cx) 充分 接近 的 曲线 y(*) 上 的 信 均 不 
大 于 QC[yotx)]， 亦 即 泛 函 的 增 晤 | 


AQ -Q[yG1-QUyCodx 0 
时 Que) E y Cr) e B RE PA 
办 区 域 妃 内》， 就 说 泛 画 Qu | 
(üAOLBGRECKIN.X Ti e X 
QUyCc)Me iil £& yr). ERUR JP 
信 ， 有 类 似 定义 。 

为 了 后 面前 需要 ， 我 们 还 得 把 “曲线 凡人 7 的 邻 域 ” 这 一 梳 
念 介绍 一 下 。 我 们 把 适 侣 不 等 式 

iux) - gr) ze 

的 一 切 蛆 线 y(2) 的 集合 只 做 曲线 gafz) 的 = 一 邻 域 @@, 在 此 我 们 还 
加 上 曲 号 yixi 及 yol 51 均 连续 上 且 均 位 于 有 界 区 域内 的 要 求 。 这 是 
“点 的 邻 域 ”这 一 概念 的 推广 。 于 是 泛 函 极 值 定义 可 鬼 述 如 下 ， 

车 在 曲线 y,(z) 的 菜 一 个 - 邻 域 中 任 肥 一 曲线 ytx)， 总 有 

4Q - QUyCGr))— QLy( rE 

WRA G EBR COOL EBOBRECK DE. 360b, AZM 
Cyr) HE yar) EGRE I LIT RE X. 

Gab EED h EE EAA RRR RRMA E 
为 基础 的 另 一 种 泛 画 极 值 定义 。 当 阳线 y(zr) 属 于 曲线 gifz) 的 这 
种 邻 域 时 ， 不 但 zz) 和 bz 的 值 很 接近 ， 而 且 它们 的 导数 傅 也 
很 接近 。 具 体 说 来 是 这 样 ， 首 先 设 yolz) 是 也 含 在 上 8 域 中 的 C!' 类 
di£gO, JL B UE 


THRE. GEB e—a M Fa 
DHr wA ipM EM, yale C) R df. 


3 


E 1 
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1) iy(r) ~ yr) e 

2) y= y «e 

3) yCc)t dé P ER BC 2E R 
Bod Reus rr. mida EIC gg 1 Bre- 2535. 

类 此 ， 可 定 尽 曲线 yotz) 的 六 Wre- S53, 

而 yotz) 讨 0 阶 e- 分 域名 上 区 中 首先 说 到 的 那 种 邻 域 。 

Xi 浮 [Ly(z)] 来 说 ， 如 果 Q[EyoltX)]j 是 QCyCr)j 在 yo(z) 
BE 1Ne-SpXEBagiAdü, ME AOC Ey cep 58 
RA. ATERACO Ey ONSE MESI ME. 

HEERES S8 BB SR BRE” . ETA PRIR IE 
FAHRE MARTELAAR “PRE” WAA WE”. 
cS E T 5n Bt rr IC mE E R EER. 

说 了 这 么 多 的 定义 ,应 当 让 读音 轻松 一 下 了 。 还 是 举 个 出 于 。 
IE, 


£i FEROUN =f vi ty r)dr 


EWE = 0, 700-1 EAP ESF AIL UR XS p, R 
BE Q BUR AM db AR vs C. | 
熟悉 数学 分 入 的 读者 一 户 即 知 QEyC2)]=| Vi y" Gode 


Thyn KE. TEMAER 上 是 求 过 点 A4CO.Q0030 B 
G,i»Bg C AE dh RPR BA OEE A DOWN A 
H ERE. yy xoBydeikxH£uBoisiBxum, 

WLgPÉRYD ER. EER AAR JLA R a ix 
PIRI, AAA BEREHER—-DOEGEQUyO 未必 
ARRALAR. WA, WREE JRERAEUL TRE 
节 中 介绍 。 


$3 mf TB ID 
光 拉 方程 是 变 分 法 中 的 基本 方程 。 在 本 字 中 我 们 对 形 怒 
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a jå = 


Qa f F(r.uCx).uy'CGr)de 
(CQ QUEUE E D de epi 


]. HH 
METE inRuEEOV-—f(xka.blbaGpeAm. X 
[/conGods— o (9) 


AHEM AARI FERRARO ) 成 立 ， 这 些 性 质 是 ; 
1》3(x3? 在 Ca ,5 上 和 有 爱 续 导数 
25 "(u)— 0—n( bh 
32 I"(x)i «eCezkfE 3& £8 3E f IE 
数 )。 
BEA. mWRIOGOfECaG.bJERNST, 
证 明 KAE, bh f uu yo 
0. SLATSUTE. IEBENGEGEJAEPES la 
因为 1(z) 在 zs 点 连续 ， 所 以 有 正 数 6 在 d dr a 
在 ， 当 |z 一 Xo 之 6 Bj, f(r)D 0 成 立 。 
现在 作 通 数 %(r)， 便 它 满足 a 
1) Gr € (G—óÓ,-,60) Hj, Pr} 0; 
2) Hr ela, ra — åkrer +ó, bint, 90r) 0; 
3) v(xMECa. 5] 上 有 连续 导数 。 
(Cr) 的 图 象 如 图 2 所 示 。 观 此 图 可 以 想见 这 种 函数 $Cx) 是 可 以 
作出 的 ， 而 且 是 很 多 的 (关于 这 一 点 网 详细 叙述 见 下 一 节 ) 。 叉 
B= APh BRICHETA PE NE, HAEN 
e> 0 nE, HRHEGACUGIEHOÉRNOA, RENG <e 。 对 于 这 
AYET), Mar 


[redre F fenar o 


这 与 4 9 ) 式 发 生 耶 盾 。 所 以 , 当 工 E(3 ,537z) 不 大 于 零 。 同 
a LERH, Arela, bh frd T3. I Wocrce(Qa,b) 
Bj, FDE, AFE r= a 5 处 都 连续 ， 所 以 Fa)=0 


4 
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—fCb), TE, Mrc[a,b], 
HEDE 


证 毕 。 

2. ”预备 定理 的 精确 论证 

细心 的 读者 会 发 现 ， 上 面 的 论 主 还 有 不 严格 之 处 。 那 就 是 函 
数 BC7) 的 作法 。 您 说 从 图 象 可 以 想见 这 种 #87) 是 很 多 和 的 ， 我 就 
不 信 ， 请 您 拿 出 一 个 满足 那么 多 条 件 的 Cx) 给 我 者 一 看 ! 这 种 说 
法 不 是 故意 抬杠 ， 而 是 表现 了 一 种 科学 精神 。 在 本 节 中 ， 我 们 来 
答复 这 个 问题 。 本 来 ， 我 们 只 发 作出 一 个 满足 上 一 节 中 所 说 的 三 
^P ZETEBU EE C m RET. 一 般 的 变 分 法 教 本 美 不 多 前 是 这 样 
BKL BÆ, EERDE EREA S RRA) 来 ! 这 
种 (zx) 不 仅 满 足 上 述 三 个 条 件 , 而 且 还 有 任意 阶 连 续 导 数 。 它 在 
现代 数学 舞 合 上 扮演 着 一 个 颇 重 要 的 人 朋 色 。 如 果 您 今后 攻读 “三 
LEA”. “有限 元 理论 ”这 些 正 在 医 厚 发 展 的 现代 理论 ， 您 经 
常会 碰见 一 个 落 数 ， 它 的 表达 式 荐 


íiz—r,2 —G 


1 rE lr, = rot) 


[ 
$a * 
i0 | :c(a,rx—8)mrC(r,Ó480,b) 


XE Pe de ET COGERET kT” (Mollifer), F8 x 
BBEG SDHEMSXTER.SO EARTE ARNE, ME 
就 来 证 明 它 满足 上 节 中 提 到 的 三 个 条 件 。 条 件 1》 和 条 人 忻 2) 明 显 
成 立 。 只 须 验 证 条 件 3)》 上 成立 ， 即 证 明 此 utz) 在 [a ,2 上 有 连续 
ERN. PARAE IRURE Cao TES t óH Ob YES TT. 
HARRY C Er + dabi Er — O MbEUSZR WEZ HR) o 
于 是 ， 问 题 归结 

8) 求 出 导数 办 ro 十 6)3 

b) iES), DIO ut (a) Gr 82, 
现在 来 求 % Cr 0) 

Harr t ij, $(x)—0 
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. lim Br) r+0) y 
二) 
Bl$Cr)fEAr.-OBBRBJES E 
pi Cret) — O0 
ToRZESÉERá, Crt å) 


lin  $(:2—-9(,-0). lim ££ j| o 
£—exotà- r—(x,Tà) pq prd z—rQY—óÓ 
lm 6 NEED ENS MM, 

t—- 00 TIS P X (x—2r,7— 3 
hm | ee . o L^ Hospitali& Jul 
EL t aa UN (在 此 用 了 pitalg£W 
2 t i 
Eme I = tire 1 
5m 4 2061 yo i 
eya 
g 21 
lim e*/*- lim -r= lim -F 
t=- QO P D RR 


SC lim (C-2eU Yo:+ 工 = 0 


o) wir. 0) 0 
于 是 p Crit = D 
还 要 证 中 
lm 'tír5—0 


Tt 


先 求 厂 极 限 ,因为 当 工 半 Xo 二 6 时 (7) 三 0 ,所 以 当 X >r tå H, 
g/(x)jm0, TË 
Tim aee 0 


了 ~ 和 双人 站 


再 求 堪 极限 。 
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a 17 。 
X eLr täi 
— — —— 
TEPENE Tk 
1 
lim  '(z)-— lime crn) sa 
+ 一 站 x 
TETEN °° 
oyat 1 ) 
E c MN (maU 
FE Cx A 
^ 
= lim ~ 2i*e'! y ó* Im 0 


p=- 0I 


Hy CTE Eat 4 处 的 去 让 被 限 都 是 霍 。 于 十 
En ye eres —Vryrà) 


PUES mS 同 理 可 知 (2) 在 xo。 一 5 处 环 连 
i. WECH, 

PEER EL bLA- BERGA, MAA E H 
Mdh “POEL ,8 上 具有 一 附 光 滑 程度 "。 现 在 我 们 
还 要 指出 一 点 ; 合用 刚才 的 论证 方法 ,可 以 证 明 名 9 Er Ed 
这 两 点 连续 ,而 (x) 在 区 间 [2,b2 上 其 余 感 显然 连 i TE A 
Wyl EL LLAREN SE SEE, RARE AAEL, b] 
上 有 “无 穷 阶 ” 的 光滑 程度 ， 这 是 把 加 (x) 称 为 “软化 子 ” 的 原 
加 之 一 《更 重要 的 原因 这 里 不 讲 了 ) 。 

我 们 现在 不 但 填补 了 预备 定理 的 证 明 中 的 漏洞 ， 而 且 由 于 
IORA, MURARE 得 到 了 比 IE SE s H 
结果 。 说 清楚 些 ， 我 们 得 到 了 : 

推广 的 预备 定理” 如果 函 数 # =f(X) 在 [ae,581 上 连 线 ， 又 
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NIS Berdar = Ü 


XH EI Hof F PEIR RI An TO BE A, IPEPE RA. 

1) "(xoMELa,.b) EH n PESA Cn AE EREHE N 
整数 ) Q0 

2) n(a)= 0 =n(6) 

35 PO «Ce A EEE N TRO 
AL. MAy = 二 ft7) 在 T9,581 上 恒 为 零 。 

请 注意 ， 在 这 里 若 呈 =1， 就 是 原来 的 预备 定理 。 

3. 最 简单 泛 廿 的 尤 拉 方 程 

现在 转 入 正题 。 我 们 来 研究 形 如 


Qty - T Fes ut rut sid 


BEA. ZEUG HR H € ERRE, 
说 得 更 清楚 一 些 就 是 ， 在 县 有 下 列 性 质 的 3X OB DIE 
Xu Xu. CE EAO Ey CoE Hg H& (&, ix se fE 
是 ， 

1) y(r)€Cfa,b 19, 

25) yCa):-y, YB my, Guy, EIEE ERAO 

25 ORRETARA BP, 

porous 2 SETPORCT BECTCHR IER 2 Dg 理 吧 。 此 定理 指 
B. ATA fOrMELaLbOB MP xr. ABHUUE (4. HLÉ" Cr E 
在 ， Wi f" 《za) 一 0。 这 是 画 数 在 内 点 存在 极 值 的 必要 条 人 镍 〔 当 导数 
fric) 。 这 种 便 导 数 为 等 昨 点 叫做 羡 数 的 驻 点 。 

世 许 您 谍 到 这 时 时 全 提 出 一 个 问题 ，“ 醋 然 泛 了 范 是 函数 法 念 
的 推广 ， 泛 函 的 变 分 挟 函 数 的 微分 的 推广 ， 而 泛 函 极 值 又 下 函数 
RERET, MA., HERAA Ea AR E” ? 如 果 


D 我们 规定 中 xz BR EEEEITUDX 3, 
© düimiECa.bo EUH 1E EnSuiTUE m. ifEC'te.b). C'a, b) 
Huib. 4A- 0,C7Qa.b.BDCIla.5., 
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您 这 样 查 的话， 我 敢 说 你 在 数 和 党 上 市 租 当 强 的 禄 象 力 !1 但 是 我 还 
要 反问 一 句 ， “您 打算 怎 胖 推 矿 号 了 本 定理 呢 2? 

“ETAR p 0 £C e A Sn b CS flc ge TA ex 53. SEE a RE 
TERT Xj. MAZA y OJE m Ry o ERIR E, 
QC 9 Cx» Eg.OoOBguxe GRETO 也 应 该 是 过 了 ”。 

HAE “PEME, -EEA i 由 这 个 绥 法 就 可 以 导 昌 都 
苦 名 的 尤 拉 方程 来 。 下 型 我 们 米 详 细 地 谈 这 个 问题 。 

定理 Flr, Y, y) 是 三 个 变量 的 连续 国 数 ， 皇 当 点 
《z 在 平面 上 的 某 个 有 界 正 下 内 , 击 久 取 任 何 值 了 时 ,PCz 3 
及 其 直到 二 阶 的 偏 导 数 《 指 计 恋 旺 :,y 及 y ps EAD 均 连 续 。 
若 在 满足 

192 y(ecrpa,b5j 

2) y(a) — pw v(5—g 

3» yO BANTTA SEES B ANA 数 业 合 中 ， 
TAC EERE ur) ER, HE BB SX 
yir) Æla, bJ ZEE SUE, PA a u CO E a E 

F,- F,=0 
du 


GEE, BrUBQUuOGO Eg GO. PREM RES TRA, E Ee 
Qly EVOBI SET 1p e- Aphk ip Boc Adis MB. EUR 
第 二 章 和 2 中 第 5 D. 

证 明 RARAC Tcl £Eu( E BUR fi Br ED vC) 
JE yx) 的 某 个 1 阶 e- gp AR, WE OLv.Cx)gz CHE 
QUu(z)), Br *vG)iyGOB1Bge-SbpP, fuori 
是 

Iu tT)— yCzx)] xe 
PHCSLEIMCNE Carela, bD 
yil )ECEa, D] 
现在 ， 任 取 一 个 函数 ?5Czh， d 
1) "(zr)€C'[a,6] 
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25 "(a)—0-n(b) 
于 是 ， 当 le Zape Wéini ulet eG) TEgCXO BU 1 
Ere- ARER Z Pg. 
事实 上 上， 于 jal 充分 小 时 ， 
yr yr = jer) Ee 
jui Cr) y Cr = Jan' a lze 
JE. yN E 
TEPRO onta ERREUR E, HILA 
QUyGr) + antz)]— yla) 
当 &z 二 0 时 
$0)= Q^ uCx)J 
Biki Ei, Mayla )/4ce — 0 处 取 极 信 。 du EL EREN 
Kn p CO DFE, Uu) 
$'C0oJ):- gD 
但 是 


P= 一 全 是 -| Fr yr)tan(e), y’ cs) 
amj e pos 


ten'GoyMs] 


Dos dOrntemeF, (Cy +a’ )] zz 

ED CF nz) F, n e)a 

=| Fcodzs Furna) | {ns)- ECFy' ydz 
=S F,- -EF yilda (nla)= 0= 7(6)) 


EP oE. FJÉv'COoDo— 0, B 


D — RRE LEE E He EIER A H E tei REB] TAR 出 x=0 
EARP (OBTER. HE, Ee Lb D SEGA— BE. GHp PE gm) + Yc 
BrMEQ—Bpe-SBhkklcPJ. E Sabte Rh Pin GEN, BIB, Ii oo EROACFPOO 
ej, piyi BHrbla- 是 定义 三 的 寺 点 ， 
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. 2j ^ 


fr —-d F'n rdr=0 
AE eI > 


REBER, (oA SESTA, BLU 
CF) E R eyy =F t Fy + Erry 


TR- F RET RES. TEF, -LF o RGEHRIRNG d 
nz) ih E AR 
1) n(r)MELs,b) EH XE SS 
2) "n(a)e 0 =7(b) 
的 任何 函数 。 由 本 章 3 中 的 预备 定理 由 可知 
F,-—— FpD (10) 


这 就 是 著名 前 尤 拉 {Euler) 方 程 。 证 毕 。 

这 个 定理 中 有 一 个 条 性 显得 不 自然 ， 就 是 隙 求 活 函 名 [CCz]] 
在 具有 二 阶 连续 导数 的 曲线 y(z) 上 取 极 值 。 我 们 说 QCylz)] 在 
yt+) 上 取 极 值 ， 是 指 取 弦 极 值 。 了 岂 就 是 和 dh 线 y(x) 的 某 个 工 阶 
5 一 邻 域 中 的 ' 类 曲线 相 比 较 ， 泛 函 久 在 ytx) 上 取 极 值 。 您 想 一 
下 ， 在 y(x5) 的 分 域 中 的 那些 曲线 都 只 看 求 是 C1! 类 的 ， 独 Hur) 
小 身 却 要 求 是 C 类 的 ， 这 就 显得 不 协调 ! 重要 的 数学 定理 往往 具 
有 内 在 的 和 谐 性 。 许 多 人 记得 大 诗人 杜甫 的 警句 ，“ 二 名 三 年 
得 ， 一 哈 双 泪 流 ”| 但 人 们 往往 不 知道 数学 家 不 把 数学 定理 琢 姜 
到 十 分 完美 的 程度 是 不 会 喷 体 的。 后 来 ， 人 终于 证 明了 :， HTA 
Qry(z)]= [ Fes y) y COM, H 

1) Fir, y y Or yky z^ ERI B fa RES, 

2) F,r'$0, 

350 Qlyir ECA hyo LRR, 
虽 吡 y(x) 有 二 阶 连 妃 导 数 。 这 个 定理 证 明 信 格 。 但 是 由 此 可 知 ， 


© 这 里 7tx) 所 满足 的 条 件 ， 比 人 3 的 预备 定理 中 的 条 忻 少 一 条 。 对 ?了 ts) ER 
BRA AT On 就 越 刻 ， 预 硅 定 理 的 结论 更 能成 立 。 
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REF sS 0， 可 以 下 消 定理 中 w(x) 是 C? 类 曲线 这 一 条 息 设 ， 
这 就 自然 得 多 了 ! 
我 们 把 满足 尤 拉 方 程 (10) 的 函数 叫做 “ 泛 AO m e” 
“ 极 带 ”这 一 稻 念 与 分 析 中 的 “ 驻 点 ”各 当 。 对 于 OX eo (Tg 
S. We'GQ)— 0 HELA OREHE E EN” ds Rer) 
TERR; HB. XE ARP” po LAQ ya NERE 
WR. iE. fk—cE/R4PT (GARTNER 8 
Qty(y2e | Flay y TEL LUCAS EZ 极 MÀ "RO 
是 极 带 。 
现在 来 淮 一 些 例子 。 
BI REALY) jt -ory +y dxi 
解 ， Fir, y, y ery ty 
S.P .— 0 
PF, -—-—2rrÓ2g! 
-d.Fy'--2:e2£* 
PE, GbR Go HI 
—2425*—0 
Br ELTE EON 
gx) 2*2 C. C. (人 | 及 CC, 为 任意 常数 ) 


例 2 求 泛 函 Q5y(z)]=| yty’ 7)dz 的 满 是 yl 0 ) 一 0， 


y(1) 一 二 的 极 带 。 


M. Eir, yy Saty- y"! 
NN o 


D xd rms gem om S WE RT. xU RODA 
BUYERS ALEISEL-IENURON. DDDTpA CUIR HD SBP? —i8, SARER 
ARRE 高 等 数学 教程 S EE EE 1:23 3 
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s UD 
"E E X 
Erud dae O 
TA. GBRXÀXDEUA 
r*—2y'-— 0 


y()o d rt EC nC, 
Xuy(0)2C,;—0 


1 1 
=> +C, +C, = -— 
y) 24 1 一世 


ed. C= 泛 函 的 极 带 为 


1 ， 了 
=i 
Hea Taa 


8i RirQUu G2 [ vides iet 


E: 人 


SAPFy'-0 
TERIERA 
A 
但 
F ,=2y 
6U—O 


WEEER y = 0 (RERO =f ydeu) 
一 上 上 取 最 小 值 ) 。 

RATER, Er, y y PERSO WE 
化 为 

Fr,y)-0 

2X, GEBÁCRIMNOOÓHRT GRWOBAGS REDEE 
程 ”， 这 句 话 是 有 毛病 的 ) ， 所 以 它 的 解 不 含 任意 常数 。 在 这 种 
情形 于 ， 龙 和 镁 方程 的 解 往 往 是 一 条 或 几 条 过 线 ， 央 此 不 硼 再 要 求 
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EF » 
它 满足 任意 给 定 的 边界 条 件 
y(a)-C, 及 2 一 
ds 
班 在 ， 我 们 再 来 研究 在 本 章 2p Ai 到 的 活力 只 fc 让 = 

| Tryde (还 要 求 y( 06) 一 0 ,801)=1) 的 极 值 问题 吧 。 

当时 借 须 于 几何 直观 在 出 总 [yz 在 过 440,07 和 Btt,17 两 点 

的 直线 上 取景 小 值 。 现 在 ， 我 们 用 克拉 的 办 法 ， 看 团结 果 如 何 。 
SELEY, y DSV ty? 
nP 


7L $57; f2 [523 


E (Zi) = i 


ny AC (C 为 任意 常数 ) 
工 十 二 


由 此 式 可 见 Czz 1 


即 y =D (DURIEGE SEG BERCAR iy nH). 

nyl) =Dr+ E 

Xy(00— 5-0 

yüu)-Dr-TE-—-1 
SQOCsCeODBIHIE RE uC(0)—0, 5C 120 1 的 极 带 是 直线 y=。 
现在 用 龙 控 方程 来 解决 变 分 党 中 上 明 车 名 的 问题 , 即 约 坦 * 见 
努 里 在 1696 年 提出 的 “ 妃 速 隆 线 问题 ”。 

最 可 降 线 问题 的 提 法 是 这 样 ， 设 平面 三 与 地 平面 重 直 ， 248 
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是 此 平面 .上任 取 的 两 点 ，4 点 的 位 置 高 于 838 点。 质点 训 在 草 力 


作用 下 沿 曲 线 4 户 由 4 点 降落 到 如 点 ， 所 
和 锯 时 间 汶 了 了 。 同 .AB 是 人 村 么 曲线 时 ， 了 时间 
了 的 值 最 小 设 质 点 在 A 点 处 的 初速 为 
E)? KERIM, ,应 当然 不 取 在 B 
点 的 正 上 方 ， 因 为 那样 门 题 就 太 简 单 了 。 
解法 ， 取 举 标 系 如 图 3， 江 质点 的 质 
Bm, REAY, XD f. W 


p e 
mu mgg "m 
n V= yy 
ES A E E r 
dS =V] ky dr 
dS — 
m p E 
ndred$ ty dr 
` vV 2gy 
TY gol (VEY dr 
‘4 2gy v 2g Y y 


zd oe A 
myy) 
由 此 可 见 ， 只 须 求 出 函数 y =y), EER 
Qtyla) z= f Vita 
FERRY ERRUER. AES HARI E 


1 
VEC, y, y'= Uu 


Us 1 -3 
《二 《一 五 ) y * 
F 


F À 二 e 
” Yyücry?) 
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DE 
Ee 

ly QE 3-1 i] 

2 V y(1- y?) 


HETRE BROKE, 585) Bd. 
Ws F-F(y0. BÀ 
F p= Fly) 


d mpo-F. AF a 


Ts 
T RJ HO FORERO 
F,—F,'y' -Fy,riy'—0 


现在 可 以 证 明 
F-y'F,.—C (CHEERED 


事实 上 
d F FE = E , EF i rF 
pra to—NH OD yt xU T sy 一 了 u” 
—u'CF ng t Frog") 
=y PF, Fou — FrinyJ]=0 


E Fy PF, -C 


将 py 代入 此 式 ， 得 


ny +y )-D ( 忆 亦 任意 常数 )》 
令 y'-lg0, M 


D 
Hm 1 +tgetp 7 


3g -Pu + cog828) 


x EE d | c Dsin2ndn "n. 2Dsin cos 9d 0. 
oyo tg im 


= gs ] T COS27 MITES 


| bee 
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: x= — Do — P sin2g4 E CE Jr ftx SD 


ii D 
| Uo i teop) 


这 就 是 克拉 方程 的 解 。 
4Z26-—m-m-—«v, WipA E 73b» 
| z-D(p-sno)- D4E 


| 
| y Li 1-—cosq) 


MMp-0BL Hr-—-o-ru. TE 
r-r(q-sing) 


y—r(i-ocse) (r=) 


mumm&o-| V ics acit eit COLE AO, WT 
à y 


ELUEBHERPANO fEL 4. B nde LARRE XX— A 
我 们 不 讲 了 。 

经 过 这 一 番 探 索 ， 我 们 才 知 道 过 4、B 两 点 的 “最 速 呼 线 ” 
T& de: — SR ie $6 2X 

让 本 节 揭 最 后 ， 我 们 来 分 绍 更 使 于 运算 的 变 分 记号 。 

RKE QU uy OJE, uNEuCXO fr tE 
óy (REX 828925C AO 式 和 第 《7)》 式 ) 。 于 是 

T(y,ày) 28Q 

X. óQ-gp'(o) CI Ca 
在 此 | ve(Qa)— QUy t aóg] 


e'() - l'CF,oy* FF, (dy)’ Jda Cip. 20) 
Tu (Sy) =y O 


Q äp = Fr (Su y yD s Ay 


He zh 46- 42. http://vib.hit.edu.cn/vibbbs 
n] 
óQ- af Fir,y, y dz 
b 
- CF iiy Foy dz G1) 


请 仔细 观察 OD 式 。 它 说 明 对 于 形 如 | FOz v y dr if GEI 


QQ@， 我 们 计算 变 分 65Q， 就 像 计算 函数 的 微分 一 样品。 
TE, 4450-0, 得 


f. Fyt FY dz o 
x 全 Foragiaz= Fdly) 
= sy Full -f ou (AFv jiz 
NT 
- t F,- -车 Pr)iy= 0 
但 óy 是 可 在 相当 广泛 的 函数 类 中 任 取 的 画 数 〈 详 情 见 本 和 节 上 


3D , BIA 


d 
Ei 
这 样 来 推导 尤 拉 方 程 ， 书 写 方式 简 清 多 了 。 今 后 我 们 将 用 这 种 办 
法 导出 各 种 泛 函 的 尤 拉 方程 。 
还 要 说 明 一 点 ， 在 这 种 推导 过 程 中 ， 我 们 用 到 了 


d(óy)-á(dy) (12) 
即 微分 符号 d 与 变 分 符号 6 可 以 互 换 。 事 实 工 
d(óg)—2(8y)'dz—-(àg')dz-ó(y'dz)-ó(dy) 
KFU (12》 式 成 立 。 这 种 二 换 法 则 给 变 分 计算 带 来 很 大 的 便利 。 


D XTiREPFG,yGD,9 (9)) sa 是 骆 变 星 , 当 f 周密 而 2 变化 则 d= Fdy 
+ Eo dy’ 
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$4 xc FG y y^ yte y Yd 型 污 
ERE JC Dr 73 RR 


f Eb EEUU GO S | F(z,y,g)dz 的 泛 函 扒 


导 了 尤 拉 方程 。 我 们 详 述 了 有 关 条 件 ， 作 了 相当 严 Me 的 论 证 。 
现在， 我 们 将 采用 上 节 最 后 所 说 的 方法 来 导出 泛 两 Q Cy = 
dz 的 尤 拉 方 程 。 对 于 各 种 条 性 只 作 太 上 略 
叙述 ， 读 者 可 以 自己 补足 这 些 条 件 。 只 要 慷 得 了 前 面世 节 的 内 
容 ， 要 详细 论证 后 文中 的 结果 是 并 不 困难 的 。 因 为 主要 的 难点 妃 
经 突破 ， 下 面 这 儿 节 都 是 推广 黑 了 。 

,首先 还 是 针对 较 近 的 目标 吧 ! 我 们 来 推导 形 如 QC y zx Y) 
- Fee yyl ydr itm RI. 

HRAC, y y DCAD, uifgmQTI AE 


Beg [ei 9 p EA 
àQ-o0 
但 是 
3Q=[ (Fayt+ Fry’ + F ráy')dz (与 (11) 式 
对 比 》 
在 上 节 中 已 指出 


b n d 
l'Foss'dz-- f ÂY dz vd 
x 


D Fiz yy UnA y y Ra AA, 3E D CE 
Tr gg D PR VEDUT d m RR E F7 
Q Py EMEUR =Y (D ya, yib = y, y'GDo = 


y LEE BOEIBEDEICE p GR AC, 
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[rutas A 


—y' Fo 


b h d 
一 Ày’ "dz CF ndr 
[eed nds 
(注意 6y 0)= 0 —d4y'CÓ)) 
i H 
本 


dco E 


dau 


L——Huyt- 


h gdz 
Far wat dz 


h d' 
-上 ÖY -Ja l dr 


^ Q= ftr. 一 -L Fo TaF.jea- 0 
和 意 3 中 所 说 的 稍 形 类 似 ， 因 为 5y 可 EAS EB GNE IEX* 
Gu S FQ Frot .F e RINHEARROR Ol) 
取 定 一 我 们 设 QIy(z73 在 zy(z) 上 取 极 仿 )、 所 以 


H^ Vu 
F,- F, Pera „>= 0® (13) 


RREZ m QCy(z)] = (zy,y detto RIS 


CORE RO) EEGRE, BEC, yy yO RAR y E 
分 光滑 ， 则 #(x) 潢 足 方程 0900 . 

一 般 谋 来 ， 用 这 种 论证 方法 可 以 证 明 ， 对 于 泛 太 QCy{T)] 二 
fF yy ® jdr, B Fx... y". y mx 
yy! .Vy 而 言 有 充分 高 阶 的 连续 偏 导 数 《 或 说 厂 充 分 光 
^ O EERE SIHH RI TÓ NS". 


振动 论坛 http://vib.hit.edu.cn/vibbbs 


. 21 * 
HO. FEROE) ERRED, Wy Coo gi m FA € 


d p d* 
n e a E 
dt j ; 
=" + 1} jt -9 (135 


It 2i £g X. S8 JE du — iB $23r f£. [dH 
8Q- ['cF y+ F py v Pysy" 


tet F, gy dx 
再 作 若 于 次 分 部 积分 ， 并 利用 
óy(a)— 0 =y la) s dy TO Ca) 
áy(5b)2 0 —8y'(b)z--2àygU Cb) 
即 可 证 明 (130 式 。 请 读者 自己 作出 证 明 。 


85 APE FU yi), usnm, 
Yx), yi(z2, yC), 7, yi Gro ]dx 
型 证 国 的 尤 拉 方 程 


现在 ， 我 们 的 你 述 可 以 和 更 简练 一 些 了 ， 这 样 读者 反而 容易 理 
解 证 明 的 指导 思想 。 


*SUCTJÉ 3d Y Cr), wo» Flr gy Ya dl, 


ydr kp BEz, Y YEY RY yit Yis 
Ui.Usd EEJEAIOETH. EET BEOQTEBUEDBy CT) yar EH 
ARE GHOEDy i ,yo 在 满足 给 定 边 值 耻 充分 光滑 的 函数 组 的 集 
A RBRSCBHR(D . My Cao 75027; F2: 


(D fROLCHu GO 4E S Een (t 


jf g'*iey se, Rat, lmp] 


^ 


n'h) se, (y " Bea 
R36 d EUR EE ILE rb s E ORC, 


FF ao 
证 ， 一 0 X. 
8Qs V (F, By e FL Sy FL òy, 


*Fu6gl t Fdyt tt Fuoyt jde 
m f Faestdzcóg Fa, -f ov (£F )ar 
oe (eene 
eQ- f (Fn E Fa)ow C. 
PINE (FL 
-—-F „jóv. Hz o 


áy,t |—1,2,:,mn ) 梧 在 相当 广泛 的 函数 类 中 性 意 选 取 。 
可 取 6y,= 0 =y =n =y A 


KETTE 
Boy ffi RF, -F WER, M 
F, -—— 
其 余 各 方程 的 证 明 辐 比 。 


$6 £L EZ IBERIA 


1. LARRI H 
W IE] PT DES QUU Cr) JE fU T ueo 1B 函 。 现 在 
3C UT S5 FORT 5 56 98 Xr B9 i. 


| mr LT 
iex 
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+ 33 * 
我 们 来 看 更 积 分 
f Z(r,y)rdy-QU7Zfír,y)] 


n 


-larai 
在 此 ， 积 分 区 域 已 是 由 不 等 式 “1; 0, PEREK E, 


Z(z,g) 是 定义 于 万 上 的 二 元 连续 表 芍 。 把 在 刀 上 连续 的 二 光 机 
数 的 集合 记 作 CLD)。 任 到 z(x,y) < CCD)，Q 的 值 唯一 确定 ， 
所 以 QCz(z,z)] 是 依赖 于 z(z LO BTEE B 

再 看 一 例 。 对 于 QUe(Cz, v))— Sup (z, 19, EIER 

| 

二 1 ) 上 的 有 界 函 数 的 集合 记 作 MCG?， 对 于 任 取 的 
Z(y)€ MC), QUZCz, vMEME— BUE, EAO ERAF 
LAG PEDE SEEATZR 

以 上 是 依赖 于 多 元 函数 的 泛 函 的 催 子 ， 至 子 -- 般 定义 ， 就 不 
必 说 了 。 


2， 用 F(x,y ,r(xy)， SS Sr Nenn 


程 
现在 我 们 来 研究 泛 函 


Qrz(z,0= 作 zf,y， z(r,y) 3, S Jazdy 
n 


fib, dE DIDCYmELA8PGXbpiS,HDS)uBE DOS 
f. FEEKS D EIE xESE SUE BU IKE, JELUSCREP hE 
fridfeC UD, XEK DHA RODEA E B HE emi, 
jg C CD) iii Je Vis L, y) loo plr y By Cm if 
集合 记 作 1(D)。 即 

R'CDy-iz(z,y ra ECHDY, 2(x,y) lon 


D SupZCRfRE EAR) bm, GRANEN bA EG, 
二 ”在 现代 数学 丸 献 中 ， 常 把 区 域 吕 的 边界 记 作 9D， 这 种 记号 来 源 于 拓扑 学。 
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z 2 d * 


—qr,y) 


BERCE Cr, 103 -ffrt y. oranes SE) dedy 


ERD pH ze. EU. GER e Cry) e P' (D, 
A 
óz—z(r,y)—zc(r,y) 
KERM yR. EA 
Bz lao=[ zz, y)— lr, y)] lən 
—g9(r,y)—e(r,y) 


— 60 
HEXRNROD p. uINUMGEDOsSX. Xia 
oz Oz 
oz oy™ 9 
mj 
ölz) ð ðz Oz Oz 
ð S ARS Ox dx e =A 
EE 
(注意 -9 二 一-S2 ange 7 的 变 分 》 
[i 


ag ien 


当 =f(zy) 及 62= 给 定时 ,， 活 函 Qrz+a:Gz] 显 然 是 = 的 函数 ， 
记 作 
eta)zQ[pzradz] 
HAZAJ, N ERRE MUAR ola) 在 ce 一 0 处 
kE CABRAmGSO, DÀ 
q'C0)— 


但 是 Qll, ide M Fi zi p Ddzdy 


Ezra zi z Rt, pdiptca:ópHRX, aH 9 二 0:89 Et, 
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* 3D * 
所 以 
eta) e Fu ota az Ptadp,q 
— 
AT ec) alfuraza F.86p4 F,6g)dxdy 5) 
2 


和 本 章 S 258 4 节 中 所 说 的 一 样 ， 可 以 证 明 若 泛 画 QCz(z,2)] 
=È Fisy, z, p, p dzdy Er OHE, W 


p'o) = ól 
Rt 8Q- f (cruza Fpi Fulq)dzdy (162 
这 是 一 个 很 方便 的 公式 《请 读者 将 此 式 积 X010 式 作 一 比较 ) . 
H9'(0)—-0, # 
8Q - fce ,82+F ,p+ roodrdy= 0 (17) 
但 基 


-E (F. z) = åz ee -2 taz) 


= 62 -S (Fas Fd (S2) 


a 


oF 


= ATI op 
; "E. t 
( 符号 5 与 -可 交换 ) 
aC F^, 
PBp= (FS2) — Eeg 
3 CF.) 
同 理 Ey 52) 2522 = + Fsg 
dF,) 
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于 是 


fersp+ FSodrdy 


B 


Ec (Foz) 62 ) jazrag 


- fex» tX CF Jszdray 
D 


IEA -E HCF 52) Jedy 


= | P ,ðzdy=F òd)" 


er 


= | (F pdy ~ Fudzyóz- 


dn 


cmo PFug)drdg 
Dp 
=- ffr CI uu XC) Jósdzdy 
I 


“jj( es -L F,.-— X Im jocos 0 


n 


n d e E HE 
B] SF. um 3 Fa -5y Et, H. Óz Hj XE T8 3 JP" i9 S] 


BARTER, D 


O ”在 这 里 用 了 格林 公式 站 (52 St asque ay uaz 
n 


e rn 
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LJ ST .- 


ð " 
uar —j-0c (18) 


F,- 
xau mo -[Irou si oodd vtt OSA 说 来 ， 


XCHRAREEAMTPEBR. EFEN, uA RE”. nmm 
(misuse pi Ddadufzz Co o) ERU GEWE, Xl 


z(xr.u)M 当 潜 足 方程 OS. 
例 考虑 泛 画 


Qcece x AC) (3, ) sev 


-ffe + qDdady 


Z(r.U) lan PCT y) 
pir y EH ETAR LEES ME AO TE BE 
zir yO ERR PAi y 应当 满足 什么 形式 的 方程 式 呢 ? 
HAF, y pag ptg 
sF acap, F.-—82q 
代入 万 拉 方 程 GO, f 


fos taurc 0 GEDAAMERID 


! z(z, Y)'an= Plr, y) (19) 


a e gez 


He~ V Jya 70 TEAR Riu E. 


D xUFEÉQGWG— r0: = | [F(z zs uum 人 a 
3i )dzedz.(U Enim Rb S, FEARR. ) 可 证 朋 ， ZERE- 2 


T "zo EREA., HsAEE HEF,- p 


> 二 
fg mt EES 
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HER OI. fie ERE EEDOD-DUOSD) EEH, TE 
JE D A elkbiu bET, TROD DERREDGEGUPEDXE s 
这 种 问题 称 为 “ 技 下 拉 斯 方程 的 第 一 边 值 问题 。 此 例 说 明 ， 求 证 


i o = (((CS- + (3)》]irdy 极 信 的 问题 可 转化 为 解 i d 


磅 斯 方程 的 第 一 边 仿 问题 。 反 过 来 ， 在 一 定 条 件 下 ， 解 拉 下 拉 斯 
方程 的 第 一 边 信 问题 也 可 转 作 为 解 上 述 泛 范 名 的 极 介 问题。 十 妃 
世纪 的 大 数学 家 黎 曙 《Riemann ) 曾 利用 这 种 观点 来 证 明 拉 十 拉 
斯 方程 的 第 一 边 值 问题 的 解 存 在 。 他 的 方法 很 有 尼 发 性 ， 可 是 在 
他 当年 的 论证 中 有 一 个 漏洞 候 。 后来， 另 一 位 数学 察 外 和 尔 斯 特 拉 
斯 CWeierstras， 此 人 可 说 是 十 九 世 纪 中 最 严格 的 数学 家 〉 指 出 
了 这 个 毛病 。 

再 举 一 例 。 对 于 泛 函 Q Cely) = ffe oan dady = 


f dz (5 zd zdy 又 z(z,) | =P= iR IE TE 


HEARR, y ERU ERRED )pioeigo 1 
ERCI, y) 0, HoMbzOO 0 (ERNA )Q[z(x,y)] 
时 最 小 值 0 。 
但 在 万 = 万 U6D 上 连续 日 在 刀 内 为 零 的 丽 数 ， 是 不 会 在 边 
*OoDErH&3dE* (quA. WMA ho- 3 是 不 能 成 立 的 。 这 就 是 说 ， 
此 泛 淫 在 C(D ) 集 合 上 不 能 达到 极 值 。 鲍 


QD 德国 数学 家 歼 曙 1851 年 在 他 的 博士 沦 交 中 论证 了 关于 共 形 映射 的 基本 定理 。 
在 这 篱 论文 中 ， 他 不 加 证 明 的 认为 有 函数 vtx,y) 使 区 画 Quan f S irua 
drduEdEÍÉ. i 
€ s (srcaodzdysbiEHsbsE S Ole (x,y) 在 集合 CCD) 中 
Dn 


&mh). RAITI., MIARE, Up c0, € CUD,Qüzr. 有 下 
id x, tHdpC tU B, QURREX:EIT MES Ok. 


振动 论坛 http://vib.hit.edu.cn/vibbbs 


* 49 æ 


$7 kx da m 


一 、 基 本 定义 
1. 泛 疯 ”把 具备 茶 各 性质 的 函数 的 集合 记 作 有如， 车 对 于 住 
Af CD, THERA — SE GIBCRUE ON ME. REE Q 7 ff 
wp fagyis. in Q--QUfgmQO-—QUfGO Ji A 
fO) ird de cpu ER, SQ — AEri y] 
2. 变 分 
《1 >》 上 趾 数 的 变 分 
XDT—27ü Eu x), MECHER 
y= YL) yx) 
在 计算 时 ， 符 号 与 5 可 了 互 换 ， 即 
d(óy)--3(ay) 
对 于 多 元 国 数 zfz,3#)。， 规 定 它 的 变 分 为 
z— 2(r,y)- £Qr,y) 


$= 9 。 在 计算 时 有 


o. 一 一 一 一 -一 
u Or P + 


C2) 泛 通 的 变 分 
o) RQ 7 (或 QCz(tz, 切 的 增 量 的 线性 主要 部 分 电 
变 分 ， 记 作 5Q。 | 
o) XimQCuGOJ- | Fer yl), y (zx)Jdzx 的 变 分 存在 
CEFER . W 
5Q=p'(0)=f (F ayt F, ay! dz 
(Elit, e€(a)- QLu(z) tady] g(R)ÓVfiig) 
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* dd0 » 


e» EEEn] fF nnen, S5, 
L 


2 ð 
Sz jdzdy 的 变 分 存在 OFRAR a Sp, P 


= 2), 
àQ s [cF.a. v F ôp+t FJóe)dzdy 


h 


3. 克拉 方程 
C1》 RTEA QUOI = | Fr[z,y(z),y'(z)Jdz， 其 
尤 拉 方 程 为 


r 28 - 
d ie dz F,-0 


(2) RFEA QUIE 人 FC, y), y CO, ns 
vO (zjadz， 其 克拉 方程 为 


CS) XETIEBUE QC Gs. (az = | FG ism 
VasyismUsddz, 其 尤 拉 方程 为 


p, a 


dr 5:(—0 (i — imn) 


(4) X-Fi Q(z(r,y2J SIDES 
HI 


2-]4zdy， 其 尤 拉 方程 为 
à 9 àz àz 
E EE E 


C 5 ) XDpiÓimRQObu, r—5.:.)] 
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a dl 。 
afe A s, "t 
2 VERTEN 
其 元 拉 方 程 为 


1] 


ch ox m 


(12 M-A- HBE A A BEN., S ED FRH po ELS ES 
ARAE, 1955. 


(23 H-:M:Teznzbasz,C-B-cbokzg,DBapmanuongoe Vicumcne- 
nue, Moc«na,1961, 

C82 MA PRERE € EHE p Wi EDBSS. 1956. 

C£3 R.Courant and D. Hiíbert, Methods of Mathemetical 
Physis,Interscienee Publishers, Inc., New York 1953, 

(52 APERIA: LHR PR F PURI” ARRA HRI 
1979 

C63 RA ADAMS, CEWE SH] Jo BOE ARRA H 版 Sb, 
1981, 
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SIE 变 分 法 在 力学 中 的 一 些 应 用 


TES — EORR ART IE RR O REEDE AR, E dE ck A 
Eu SeEQgeXmERAEE pERGE (BO OLARI EHE PM EE 
Bo .mpuxe xg Jodie; 2. GEB gef dick 5 — A 
kB, AÉARBHEIEORASS€E. xxppigmis. C did un Ws Jb f 
HIER UIRTGÉBSBS mom 5M EHR 
4 WE MRO EKAR LB A Q AD AKEE dE S 
小 值 ? HOESBEEBPÍBHEGEXEGE, MAI AE Dixcep]es, pid 
谈 变 分 法 的 应 用 。 但 是 ， 对 初学 者 来 说 ， 和 作者 认为 在 学 了 尤 拉力 
程 之 后 ， 还 是 先 接触 应 用 问题 为 好 。 理 由 是 ; 

1. 广大 的 工程 技术 大 员 关心 的 是 用 变 分 法 解决 实际 问题 。 
对 这 些 同志 来 说 只要 掌握 了 本 书 第 一 、 二 章 的 内 容 ， 器 可 以 四 
于 实际 了 。 

2. 从 变 分 学 的 发 展 吕 来 看 也 是 这 样 。 尤 拉 和 拉 格 兰 日 的 工 
I REBTTSHSUKIEXT GIGS Bg. È 
拉 是 杠 1736 年 发 表 他 的 成 果 的 ， 共 中 有 后 来 命名 为 龙 科 方程 的 壮 
各 方程 式 。 在 此 后 七 十 余年 间 ， 数 学 家 们 长 期 探索 泛 国 驴 取 极 值 
的 充分 条 件 ， 设 有 取得 多 少 成 绩 。1807 年 雅 可 必 (XJacobD 在 这 
方面 稍 有 进展 、 得 他 对 问题 的 叙述 是 不 准确 的 。 直 到 1879 年 ， 外 
尔 斯 特 拉 斯 才 焉 确 提 出 六 严 咎 论证 了 泛 由 取 极 值 的 充分 条 件 。 这 
已 经 荐 尤 拉 方程 出 现 后 的 第 一 吾 四 十 三 个 年头 了 1! 尽管 充分 条 件 
的 建立 是 这 祥 晚 ， 但 是 并 没有 妨 得 变 分 法 在 这 一 个 半 世 纪 中 广泛 
下 用 于 实际 。 

3. 莹 函 极 町 存在 的 充分 条 忻 的 论证 是 相当 轴 难 的 《本 书 限 
于 籍 幅 ， 不 谈 这 个 充分 条 件 了 ) 。 

本 章 主 村 是 介绍 变 分 法 在 力学 中 网 应 州 ， 其 他 方面 路 有 涉 
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=. d3 >» 
及 。 恋 者 理解 指 学 思 奶 之 扩 。 自 然 可 应 用 到 更 广泛 的 领域 中 去 。 


$1 Nix DIRE e s A RR 


1. 哈密 尔 顿 原型 

当 和 牛顿 建立 了 以 三 大 定律 吕 万 有 引力 定律 为 基础 的 力学 理论 
以 后 ， 无 煞 的 和 然 现象 邦 得 到 了 定 莉 的 说 明 ,事情 似乎 已 缀 很 完善 
了 。 人 但 科学 家 是 永 不 漠 足 于 已 有 成 就 的 人 ， 届 人们 仍然 继续 探索 宇 
密 的 奥秘 。 后 来 ， 拉 格 兰 日 在 十 八 世 纪 旬 出 一 个 变 分 原理 ， 从 这 
个 原理 出 发 ， 运 用 变 分 法 ， 能 够 十 分 方便 电解 次 力 学 问题 ， 并 且 
还 可 以 摊 导 出 力学 中 的 主要 定律 。 在 他 鞭 的 芭 分 析 力 学 》 中 和 叙述 
了 这 些 成果 。 据 说 他 讲 过 这 人 么 一句 证 ， “在 我 的 《分 析 力 学 中 
没有 一 张 图 ”。 他 还 创 痰 了 “ 拉 格 兰 日 运动 方程 ”， 它 比 和 牛顿 的 
运动 方程 适应 的 范围 更 广 ， 用 起 来 也 往往 更 方便 。 以 后 ， 哈 密 乐 
Bi (Hamilton) 多 发 展 了 拉 履 兰 日 的 理论 ， 他 在 1834 年 提出 了 了 
一 个 车 名 的 党 理 GERAN ERWEE ， 大 意 是 ; 

在 质点 《其 至 是 质点 系 或 物体 )》 的 一 苛 可 能 的 运动 中 ， 真 实 
的 运动 应 使 积分 


s-["cr- uidi 
tl 


取 极 值 。 在 此 ， 了 7 了 和 如 分 别 是 动能 和 位 能 ，#, 和 :是 两 个 任 取 的 
时 刻 。 

这 个 肾 理 竟 成 了 力学 中 的 基本 大 理 。 以 它 为 基础 ， 可 导出 牛 
B Kx inb. JAMS EEFE, DEUS x 
m. GARENA 力学 了》 一 书 CTED., AEREA 
会 感到 困难 的 ) 。 现 在 我 们 来 把 此 原理 说 得 户 消 楚 一 些 。 

哈密 尔 顿 永 惠 的 精确 表述 是 ， 假 定 在 1 二 二 及 1 二 tt 时 刻 质点 
的 位 置 已 分 曾 确 定 在 A 点 及 BB 上 点， 那么 质点 诡 动 的 真实 铺 道 及 如 
HE. BUB 
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s-[*cr-ut (SLE “EHR” O 


PUR (8, EI T RU eO Ht AR 80 BE AAE. PRA AE AA 
d E HAI T =le T r EEEE ER, t ERE R, T S 
的 曲线 表示 真实 轨道 。 这 是 示意 图 。 

A 试 由 队 密 尔 顿 原理 导出 后 上 e 
第 一 定律 

对 于 自由 质点 ， 动 能 

m E + p T 
=T (rtty z*) 


位 能 U 恒 为 零 。 于 是 作用 量 


S Zatti d 
Sii, BOROLDTA4Ca.)5:406)0)8, 
Sici WgufrTB5Bt36.5:,c04. Hl 
rli =a; ICSSE DP 
ylti) =b, ylit) = 0; 
z(i)-—c, 2C) =C, 
显然 psf Zat ytd m SD CD (02 
i 


是 轨道 


| x-—rít) 
0g m yit) 
| £&-zcii) 


BUTS. KER ES 5 中 所 说 的 Firul) yE), 
YILT), v. yat) Jdz €) it p. dig5n[Am, FAS EH 


(D ztD,yOQD,2€00 m HB E Sen Ge ya) SET IRE BU Ro Fs y, 288l 


piz, E. DL. XORAASSDARUSAR Cb ARGUS ARO ES. 
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* Ah - 
| zc) 
道 ( 9 =y) ERRE, Me) y Cz CON BILE 
(oz (i) 
d 
F zm x a = 0 
d B 
| phy ur 二 人 
| do m 
(BEEF mL C zt) 
. F d wr 
S.F.—0 Fi-—mr "up Poe Habe p 


MAr-0 同 理 得 y= 二 0，z 一 0 
可 见 自 由 质点 的 加 速度 为 堆 ， 所 以 其 速度 为 常量 《说 得 更 清楚 一 
此 ， 速 广 是 常 失 景 )》。 也 就 是 说 ， 当 质点 不 受 外 力作 用 时 ， 它 或 
租 止 ， 或 作 色 速 直线 运动 。 我 们 由 哈密 尔 顿 原理 导出 了 牛顿 第 一 
定律。 

B) 试 由 哈密 尔 顿 原理 导出 质点 在 重力 场 9 中 的 运动 方程 

在 地 球 表面 附近 ， 可 认为 重力 加 速度 9 KETE. HERA 
原点 O 放 在 她 下面 上 , 取 = 轴 竖 直 向 上 。 质 点 的 动能 人 一 于 (zz 二 扩 
£27) WU —moz, PEER CHI * (ERI ) 


s«l (ons y? = £)-mgz |t 
t] 
F7 Ga g!-i-mgz) 


Es FPF.=0=F, FF,=—myg 
Fi-mi, Pu-mg, Pi-mz 


Q dEHIREEHuIESSASZAÓE. 
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.- AG -* 
出 F,- mio fuoco 
由 F,—-Ž-Fi=0 Fy = 0 


HI F,- F:=0 4H — mg- mł= 0 ..z——9g9 


TER A (E E 7) Yn rm Pai OE 


例 — he ERRE S HL E TERR 2| A Zap "P ET Sl 2; T8. 

HIRAET OO, MARERE E AEA m 
RM, HORGIAS ERO CO), usen zt1)) ， 于 是 动能 及 位 能 
分 别 是 


7 一 本 (如 二 9^-r At) 
ITI 一 了 i4 21312155 
UCM ocu s r= rai ytz 


作用 量 为 
s- ons $^--2?)4 KP Vut 


-[ FO aC (Y eC 2G GOD (3t 


Merata KU 
dias -4$*0-3*) K- 


A mk mK 0 
On £567  Pü-me 

mA | 
F,=— ES y F,=mý 


D mE DERRAT ARAR R ET E e EE TA 
fd SA SARRE C EMEN aiee s d ARET Ae ERAT 一 
wi. AR4]màn AUG Se dn GRÉ EX, y, zlITnEI8]AVE[BIICT XD 3E EOCAE. 
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F MK 


"E 


z F:=må 


JUR O4) TEN EE BUSR S173 P a 55 277; Ee 


t=- KZ 
| / 
t ġ#=-K t 
| i 
| z= Km 


2.  WREXGRBUJI ABAU AEE 

傅 定 原点 在 平面 上 的 在 置 ， 和 需要 两 个 吾 标 ， 这 是 常识 。 但 是 
请 注 蕊 ， 这 两 个 华 标 不 一 定 是 直 儿 坐标 。 例 如， 在 平面 上 到 极 华 
WR MAME HEr (fS) Ho CREDO 所 确定 。 由 


: t -—-rcoosgg 
Uy = rane 


RIIE. PRA ASI LEER 

—üxmck, PAATE on «n. H 

T— ELG Qa) 

y= Ytg Qa) CHH Ego go ap- ECT. v» 
WAG qe D AE EER. 

对 于 空 册 中 的 质点 ， 它 的 直角 仅 标 是 《7 ，y，z》 ， 若 有 独 
Xp EEq., qi, fa H 


| see Qu Qi) 
) y= Yq Qs Ga) 
| z 2 {C01 "E qi) 


CEL q,, gO aE E r, yD 则 称 9,，9;，9s 是 
z= [Bl H E REN S AES. 

总 定 空间 中 一 个 质点 的 位 置 要 三 个 广 尽 坐标。 要 确定 产 个 质 
点 的 位 置 常 需 要 3 访 个 广义 坐标 。 一 般 说 来 ， 能 够 确定 体 系 《 质 
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点 系 或 刚体 ) 位 置 的 天 个 独立 参数 9i，9z，…9r 称 为 此 体系 的 
«y^ S AAERU . WEAN se g, di, tex 98223] XD iE 
BE. SOR IR AE RA G AE 

3. 拉 格 兰 日 方程 

现在 ， 我 们 来 于 出 描写 质点 运动 的 控 格 兰 日 方程 。 先 设 帮 点 
HASAT AER] 《这 时 质点 的 自由 度 为 1》 。 AARAA 
由 广 闵 坐标 9(+) 决 定 ， 所 以 动能 了 和 位 能 0 十 9 Ee poi. d 
T- Umui HM, in 

T —U-L-—-L(t, q, 4d) 

于 是 ， 质 点 的 作用 量 为 


S =f LG, g; ġidi 
1 


因为 3 吏 极 值 O， 所 以 真实 轨道 9(1 7) 满员 


oL d dabh 


RRENJE RAHMET. J SARAR CEAk 
MAREA BUfrTEEREUXASERS s qe |, q 决定 《这 
时 称 质 点 系 有 上 久 个 自由 度 ) Hga DR glt) G= 1, K) 
均 已 给 定 《 即 在 (1,28 (f =h 两 时刻， 体系 的 位 置 蕊 给 定 ) , 当 
质点 系 由 时刻 的 位 置 变 到 :时 刻 的 位 置 时 ， 作 用 量 


S=) LG, 0 ge dn gdt 


RRN. SERTE SHIR Fn vle) e, yka), 


gir) eegal ldr H. MEAAG g, G Ci o1, 
we A) ism 
9L mc. QE. o Ci = l,e, Xx) (21) 


KREA RREH TEH, CEE EIRFA, 


D 指 当 g Ci) = 9 gpg AEH, MUESUR EGRE. 0 R00 GER 
WE ge g gine gph SIEGE. 
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系 的 运动 方程 ， 表 过 了 质点 系 运 友 的 一 般 绕 健 。 

读者 也 许 会 问 ， 拉 格 兰 日 方程 究竟 各 ， 
些 什 么 优点 呢 ? 我 们 还 是 举例 来 说 明 问 是 
HIE 1 

例 在 长 为 了 的 细 丝 《其 质量 路 点 不 
TD) xd vA. Gi Hm. 4255 
-一端 固定 ， 置 重 方 场 中 。 写 出 质点 的 运动 
方程 。 

RERAMA 5 所 示 。 此 质点 的 位 置 显然 由 一 个 广义 坐标 € 
Wax, AMET musime, p= 0 时 ， 取 质点 的 位 能 


为 零 ， 则 在 9 处 质点 的 个 能 为 
U —mg« AC —mgl( 1 — cop) 
qtti 25 Hag $679 


peqety =L mpy — magit t —cosug) 


ra 


于 点 的 拉 格 兰 日 方程 (90) np 
ðL — d ðL 


èp df oe 
由 
ÖL z d aL d à - 
一 一 — — Imgsing Dod —£—Dpe——Um vsu 
p di ( ae ) dt 
可 得 
p eT sinp- g (22) 


这 就 是 单 提 的 运动 方程 GEARRE), RERED 
阶 常 微 分 方程 。 落 19 | 充分 小 ， 用 sinz = 一 吧 取 


p 4 p= 0 (23) 


代替 上 述 方程 来 描述 单 押 运动， 这 是 线性 二 阶 常 微分 方程 。 它 的 
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p= e COS 7 十 casin i- f 


v€(0)20  9(0)—4À 


p= Ay Lin y i: 1 


可 见 当 质点 作 微小 摆动 时 ， 其 周期 是 3x 站， 但 当 摆 幅 二 大 ， 


g 
其 周期 就 不 是 此 数 了 。 

闭 才 我 们 对 摆动 作 了 一 番 分 析 。 由 此 例 已 经 能 看 出 拉 格 兰 日 
方程 的 巨大 优越 性 。 它 适用 于 各 和 形式 的 坐标 系 。 在 解 应 用 题 时 
不 必 商 力 的 分 解 图 ， 只 须 写 册 动 能 人 及 位 能 尽 关于 广义 坐标 4, 及 
广义 速度 ?9，( = 1 ，…，) 的 表达 式 ， 然 后 代入 拉 格 兰 日 方 
程 即 可 。 

在 这 里 我 们 还 要 说 明 ， 因 为 本 书 只 是 简略 地 介绍 一 下 变 分 法 
在 力学 中 的 应 用 ， 所 以 对 于 拉 阁 兰 日 方程 成 立 的 条 件 《 即 对 什么 
样 的 质点 系 ， 它 的 运动 方程 可 用 拉 格 兰 日 方程 描写 DUBIE 
论述 。 读 者 欲 希 其 详 ， 清 参阅 盟 道 所 著 《 力 学》 或 国内 近年 来 出 
ESET 


$2. 变 分 法 在 力学 中 的 应 用 举例 


在 处 理 各 种 工程 技术 问题 时 ， 常 常 需 到 推导 出 描写 基 些 物体 
运动 的 微分 方程 ， 然 后 利用 微分 方程 的 解 来 找 出 运动 规律 。 这 个 
办 法 是 生 顿 创立 负 。 他 找到 了 描写 天 体 运 动 {后 来 看 出 还 适 凡 于 
各 种 物体 ) 的 微分 方程 ， 然 后 断定 天 体 运 动 轨道 是 圆锥 曲线 。 另 
一 个 著名 的 例子 是 麦克 士 威 (OMaxwelD 找到 了 描写 电磁 场 变 化 
的 微分 方程 ， 然 后 根据 此 方程 断定 电磁 波 存 在 ， 这 是 1868 年 的 
事 。 直 到 1887 年 (也 就 是 在 麦克 士 威 根据 油分 方程 作出 他 的 科学 
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预言 的 十 九 年 之 后 ) dhik Cleto 村 用 实验 证 实 了 电磁 波 的 TA 
在 。 在 科学 忠 上 还 有 有 大量 的 例 注 可 以 说 明 微 分 方程 的 巨 
大 作用 。 在 我 国 工程 界 ， 现 在 大 家 习惯 于 用 达 朗 酝 尔 原 
理 列 出 力 平衡 方程 ， 由 此 得 到 运动 方 和 可。 这 种 方法 园 然 
有 和 它 的 做 成 ， 但 是 着 用 变 分 原 惠 来 恒 出 方 穆 ， 却 有 更 大 
的 优越 性 ， 在 本 和 节 中 ， 我 们 将 从 各 方面 举例 来 说 明 这 种 D 
方法 。 
MEE XE C EUREN - 
RAME 6 BUR. WR Lou rika fzgbdu P" € 
u(z, 1 MEWE, Ie! &mup3gEJIF 的 作用 下 CHE HEE 


行 ) 所 发 生 的 长 度 改 EH ke EREADER, SAH 
的 横断 面 面 积 。 于 是 相对 伸 长 
zi E UM a PESE 


MERAT ERE PED T Br. 
Jul 


goXxtem u(rcdrt)y—-r—u(r,t) 


zd 
m Ait a 
Ar 
eux riri) OLALI) 
Fir- Ar, t)j—FES5u,(rct82mx, t) 
在 此 式 中 ， 令 z- 0， 得 
Fir, ff Su, tr, f) 
Fiz Ar, D= ESuy rci, 1) 


现在 求 弹 幅 位 能 已。 由 弹性 力学 可 因 ， 漳 性 体 T+A 
在 谈 形 过 程 中 克服 内 力 所 作 的 功 ， 作 为 能 量 积 莹 在 
体内 ， 此 能 量 称 为 弹性 位 能 。 村 眉 Cr, rAr] 
所 受 弹 性 力 为 ( 见 图 7》 国 了 


Ft) 
A 


Fas axt) 
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* E LJ 
SEu,QGr d dz, ty SEu {rt S Eu (Cr üzfs tir 
CO « g« 1) 
于 是 ， 在 di1 时 间 内 ， 此 为 所 作 的 功 为 
Eu, (rr tr wdt- Eu, (rx, t)u,dvdt 
| — Jd | HH 
BEIE (x, rcx) PARRES. JXCTOGRÓEBEGEBR 内 
力 所 作 的 功 ， 当 然 与 此 式 反 号。 所 以 在 Cr, f024£3 qp] dm] Hé 
PREAH E: Cr, r-ocfr 中 的 位 能 微 老 是 
dU/—-—SEu,(Qr, f)u,dxdt 
于 是 ， 弹 性 位 能 为 
U = -| | S Eu, a, udadi 
RE 


— "dud eh E 
scan KL MN Put) usns, a 


——SE INE LIT yt 


SE RA uu drdi 


— 十 i i Q t 2 
=5F KZ | Ll y 
=SE | Luidzc 
SE | dz (24) 
BbRGD8B, Ær E (x, r+Arl 中 ， 动 能 微 元 是 
dT ——pG)utiG, Sdr 《kz) 是 质量 密度 ) 
所 以 动能 是 
T =| xeGoSWiG, Hda (25) 
23 


拉 格 兰 日 函数 为 
L- r-U-|'(ocosvi 一 二 3 已 ze E 


T £ifÓESWGIEH, HEUTu.C0. D 0, n,U, 029 Eu, 0020 
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. p3 æ 
VE IR RES 
s [*cr- vai 
=[ Hesa: --4S Eu? Jordt 
Ti, SE, f)Bg its, BI 
Stutz, 07 [/ [FG t ens Geb, 
1J 0 
u.(r.t)Idzdt 
这 是 在 第 二 章 86 中 讨论 过 药 多 元 函数 的 泛 函 ， 它 的 尤 拉 方 程 ( 参 
AUH Æ 
à à m E = 
F,—-2 F,- F= 0 (u= P. 5,0) 


F-Le(G)usS -lSEwuo lo(ose -15Ep 


Ag" P. seo F,-—SEn--SEu, 
F,—p(x)Sc-—p(z)Su, 


ð : E = 
3 OE ar (0)9)-0 


uum Eye. C301 TTL KE I (26) 


当 ptz) 兰 p【 常 数 时 ， ang 


u, ET (27) 
XX E la 0 2 77 E P n n sx te sg p e 
LEE Xr qEERSDE EE 
《 设 弦 两 端点 固定 ) 
PERR., Mx ur, 
1 表示 在 时 刻 + 点 + 的 横向 位 移 。 
HOHER (x, tAr], (ERE 
过 变形 后 移 长 度 是 Bos 
di-Vi-cruidz, 所 以 此 段 的 伸 长 量 是 


ii 
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Viludr—-dre(Vicul —1yfr 
SEERE ) dr= z uidz 
2 2 


Quuicl. Walik 68 €) 
Eiet, iget R RAMIE., BI 


dti -: B uida (28) 
于 是 位 能 B 
U=| 2 T (29) 
X Si Bb 
T=| iocus COC GE BUR SERO — (0) 
fE FH REOS 


fef f 3 1 » 
s=| oe D- d Ku, P) Jardi 
RESTÉ, £0 HZA. CAA 8 CBDJSF3018)) 是 
à 2 "E à * - + -= 
Ts t) Es [oCrj)u;(mr,2)] OQ 
A 


A" HH. D: 


pr) 
这 就 是 方 田 (26) ， 当 P(x) 三 Pp( 常 数 )， Jr mat, 我 们 又 得 到 


方程 (27》， 可 见方 程 (206) 和 (127) 上 县 可 以 描写 细 杆 的 纵 振 动 ， 
X up £859 25 d p Dez 

例 io AE ESI T. 
问题 前 条 性 与 前 例 基 本 相同 。 雁 一 区 别 是 设 茂 上 各 点 均 有 外 为 作 
Hj. st) EE Wp dpfr EE O rt E B. 外 力 的 线 密 度 为 
gOr,1)0D, EAE BFITESEERU X ru rER BATES. dr, 
TERTA t, UE LAr ABBU, D CGEgedgpSEER B, 5RÀ 
作 的 动 可 近似 地 表 为 

— ptr ulr, dz 

(OX E EFH SE, Bampa t, 


rz 
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末 以 ， 在 了 时 刻 ， 外 力 位 能 为 
- vG, futr, Ddr 


这 份 能 及 也 积 昔 在 弹性 体 中 。 可 罗 防 的 总 位 能 应 当 是 弹性 位 能 与 
SD iez, Ti 

U - | CER utes D - $6, Iur, D 

T- | Lptoywtte tdt 

S 一 -|. "ue Lou -lKueyu Medi 

-Lv 1 Bagues loli LO RR 

F-.pui ;Kus pu zoa AP tu 

` Ó n ME à - E 
由 Pues od, px 


vir.0)— pGaDu, t Ku,,-— 


. mue a qo Nar. 

M U, = IES u, {1 jS (312 
K (xf 

Mpzso(xyssoCEBO RI, eat, HoD opes d 


a DOEARPCGHUSgsDqwgu, Weiz, 0) 2g EARR 

. H,,—aG^5u,, - Vix. f) (32) 
方程 Cay 及 G2) — i8 ap Jp $E. 

例 推导 如 的 自由 振动 方程 

取 坐 标 系 如 图 9。 uir, y, 
O 表示 在 1 时 刻 腊 上 点 Cz, 
y) 的 位 移 。 在 膜 上 任 取 一 小 
H, EMERY 
4S -Vituirui-1ldzdy 


~ (uicuiMrdy 


Cuimil luu S E 9 
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于 是 ， 在 此 小 块 上 弹性 位 能 是 
dU ui uiydrdy (KREYO 


Rn ye 十 udrdy 
又 T- | ee y)u?dxdy 《plzx,y) 是 腊 的 质量 面 密度 》 


作用 量 为 
S «Jat [osi -Kca Tul) juzay 


=j [FE yt) ur yt), 
£a 


u,(r,y,t)ldxdgydt 
RESA, y, NDA, ry RAZA, CHEE 
方程 〈 见 方程 (187)7 是 


Ó 
Peras Ea dy jew jac? 


S F — p, yu A Cite ub) 
*n F.—0, F, ——Ku, F, =- Ku,, F., —pu, 


á ð ð E^ = 
cg uty An Br PT )= 0 


天 
"E dues V) dn PUR) (33) 


Mpo(r, =P 50, $—-—a. TË 


u,,—a?(u,, t uuu) (40 
这 就 是 膜 振动 方程 
例 TRA H ESO E 
取 坐 标 系 如 图 10。 Hul, HE 
a: 时 刻 梁 上 的 点 * 处 钢 位 移 。 由 弹 
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æ 后 了 " 


性 力学 可 知 其 位 能 与 曲率 和 区 平方 成 正比 。 于 是 ， 在 梁 段 fx,， x+ 


dU =A -Muüce Kuzdzo 
| ] 


Hi 
(13 uii 
E OU 一 i Kut tr, idg 

J D 
X T= 二 | PLu e, dr (OC EM AARE JE) 
"n 
^C FREA 
Ss=|, | (eu — Kot rdt 
- | | FC tunis d edt 
ri u 


Fal eGyutG t) Kul (e,t) 


HTRAS, CMÈiH AAE 


__9 - : 
Is Ja amm m da Hos cpu, Pep (35) 


在 第 二 章 中 ， 没 有 讨论 过 这 种 活 函 的 尤 拉 方 程 。 现 在 对 方程 
(35) 作 一 简略 的 论证 。 
jun] Hg uH pgs — 9h. A IBIZIN 


S =|| Ftz,y, u(x,y),U uH, Wr ty ld rdg 
Ù 


在 此 ， 区 域 中 的 边界 3D 及 函数 五 均 充分 光 油 。 由 585 二 0， 可 得 
泛 函 S 的 龙 拉 方程 P 。 推 导 过 程 大 致 如 下 ， 


es |l or jugi. 8p 4-9. ôq +E ar 
D 


op oq dr 


aF 
ðs 


ås Jazdy 


Q 由 弹性 力学 可 知 ， 常数 六 = 去 五 了 ， EREEME, JE EE NEUS D 转 动 
mE. 
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TEIE, H= p, uQ,—Qg, S, —r. Hp =s 
在 第 二 章 6 中 各 指出 
中 Fap+ Fð dedy 
D 


三 [GO +2 F, Jaudzdy 


dy 
X 
Eor Eo yn 
-fje eS SEa ds de 
— HC : " (Su, Tz € (Óu, j Jdzdy 
Eo "x 5.5 M a (3 in. eso 
Aic: iE. Su, + E! t3 y Âu, Jardy 
用 格林 公式 得 到 
[Ls Eae.) GE. een 
gl STan, )ay- (Eou, jiz J= 0 
Cru Eu EUR EE, S898, ™ 0 —8u,l ao 
x JBE a P) ye 


D ERJA uana, slln-V. uilan-co,. TEREARIXPDIT[A X PRH 
大分 光 许 的 函 烙 类 中 研究 汉 鲜 S 的 极 值 。 设 SE utr, wg) ERRE, fixo 
uir, um Wm HUE. 


Ean si nanmi 
t 
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. 59 > 
E le s uod n 
人 


MG bau T 


但 
j Cs VIS 7 ( zE Ou 外: zdy 
= Ie zr" N E òu jde }=0 
nex, 


Je SP ar, 3958s yo (dy 


-JUE (GE) (Eora 


r 
ace] SE ò oF 2 3F 


人 


F Qu Or Jp ðy Og 


Q^ OF, o dF 
ort r oy? 


ES A CR XELE E E PER e EAE EE FU u ARR E. 所 


A Jaudzdu- 0 


oF ð oF | OF, ð dF ð oF 


一 一 一 -一 e o ammo UMP i niai i m. — -——- 


Qu ðr Öp oy T dx" o? Oy* Gs 


= 0 
(36) 
3 $545 =|| Cy nn sss ss sod cd tti H 


方程 《此 方程 用 得 不 多 ， 在 第 二 章 中 玉 避 到 ) 。 
BETE, HTA 
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t. 了 
s=| | Jig d Na IE NI IR... Mag 
"idu 


TR EAA tE EGD XA 4Apro- d patur ped dcr 


TE d poU p ; mE. 
KE. PS. Pr FS, (£A (OD X. fi 
Vn _ 3? 
pt | dis oh, grum u (357) 

at CQ. 


这 就 是 梁 的 百 由 CH ika Fi 
例 ” 求 在 重力 场 中 定 轩 的 轴线 的 — EILILITL—-—————x 

JER Gr AERE E HL AAT 75 on : 

程度 很 小 ) ,注意 现在 深 受 理 力 作 崔 ， 

Br ABE RAE O HIRA ADS XE 

BÉABOHTOEÓRdRÉRGR. GXXWIPGUIGEYJUE. GA , 

ASA AnEM1. Bb o s(x) ETE 

E roibit GAE u 与 时 间 无 

XO. YE Lert Az) Em RA RE = 一 Gatrdre MEER 


的 弹性 位 能 ~ 二 73dz。 Laut LAA, M 


村 11 


| lE E 
v=| (; Efui —-gpQr)u(x) jus 
HERZEN 
as LE "e | 
S- MNCLIDS geG)uGn) edt 
F-LEIutGG) gPa) 


Jb i27; EE 7E 
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" Bl] * 


sa GPJ EAn H = À 


uq =L (38) 


这 是 4 阶 党 往 分 方程 ， 它 的 和 分 曲 级 在 4 个 任意 常数 。 它们 可 南 
HARPEN. EEE fRA EJ R5. 

这 个 例子 也 可 用 着 下 方法 求解 。 妆 庚 处 在 平 险 在 (稳定 平 
8D 时 ， 动 能 为 零 ， 位 能 最 小 ， 所 以 


U = '( : FEFat,.—090(rYu(xr) E 


V DM 上 
=f Pr, uir), uiy lde 
f 


取 极 小 值 。 

XEM UCuCx)HE-— AG R dU fe uCr) EHE 
ii, Weto 满足 尤 拉 方 程 Q2. 
写 出 的 方程 仍然 和 刚才 的 相同 。 

例 设 梁 的 两 端 图 定 ， 在 梁 上 某 
点 有 一 集中 荷载 。 试 用 广义 坐标 求 浴 
HEA Y 

HA ERIEGHEBISHDTR. EREE 
IRMCBREH—EBAPPPH. xxt 3B 
An£bnH yxp Hr (ourien 2% 


. i HI 
ET cet a SI x uc "€ 


y= asin = + a, sin 7 
了 世 就 是 说 y-—Jf(a, a4 c 0, 0X 
现在 把 Ca, Gp es Gu …) 视 为 此 质点 系 《 即 梁 ) 49/7 E 
标 。 于 古 此 系 有 无 穷 个 广 只 坐标 。 过 去 我 们 只 衣 论 过 具有 有 了 眼 个 
广义 泽 标 的 质点 系 ， 但 是 现在 我 们 不 妨 继续 启 论 下 去 。 

C 点 处 的 位 移 为 


yc) a, sin dS 


3L HORIS 
B ieai in 一 


le 
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* 62 + 
外 为 产生 的 位 能 为 人 
nk e) ~ P A asin E 
弹性 位 能 是 
U,-— EJ IG " dr 
: s n - " 
TE! {È T a,sin x ds 
n* " 
i Sntas® 
h-i 


于 是 总 位 能 为 


^ EI : 
U -U,« U,— — P "S osin Lre ET N tat 


=U (a, M. Guy ae] 
为 为 鞭 处 在 稳定 平衡 态 ， 所 以 忆 段 极 小 值 ， 于 是 


ðU 


EP Ü (12z]|,2,:-) 
RA T into. Psin— = 0 
* a 2IP Cn ncc 
-e "Ca EJni i 
由 此 得 轴线 方程 


TO uwp£kx HEAR. MUEAGENSBEETSTE:5—T 


ù w yAn 
BEP 754-8: ES T sin E sin PTT dr = 和 
n I i FI 当天 站 
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P "ncc ni 


^ 328 " 
Y= DEI n7; 


由 此 方程 可 见 ， 当 一 0，(! 时 ，“ 二 0 和。 这 就 是 说 ， 当 集中 力 
ERER T EE SB, 


$ wG x à 


C12 J.J/I.JIaggay, E.M.JIudaongzm, € JE. SARAHA 
[P 1559, 
1223 BHE d EDERA Jo FEHR 1950, 
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第 四 意 ” 变 分 法 的 新 发 展 一 一 有 限 元 法 


在 第 三 章 中 我 们 举 了 许多 例子 ， 对 于 它们 我 们 都 是 写 进 有 关 
泛 庚 欧 尤 拉 方 程 就 止步 。 但 是 在 解决 实际 问题 时 ， 还 要 求 拒 这 些 
方程 解 济 米 。 令 人 坦 往 的 是 这 些微 分 力 程 多 数 是 不 能 积分 成 初等 
明 数 的 。 计 古人 们 只 好 另 想 办 法 。 所 请 “ 解 变 分 问题 的 直接 法 7 
就 是 适应 这 种 需要 而 产生 的 。 尤 增 开 始 研 究 变 分 问题 时 ， 就 用 过 
“有 限 差 分 法 ”， 这 是 “直接 法 ”中 的 一 种 重要 方法 。 记 以 无 拉 
也是 变 分 问题 中 直接 法 的 创始 人 。 司 来 你 出 现 了 里 兹 (Kitz) 法 ， 
傣 辽 金 CTanepxarn) 法 等 重要 方法 。 在 这 些 方 法 的 基础 上 ，?” E 
了 “有 限 元 法 ”。 这 种 方法 最 时 见于 柯 关 (Courant)1943 生 的 论 
文 。 他 在 这 篇 文章 中 ， 首 次 用 了 一 组 “三 角形 元 素 ” 来 解决 力学 
中 的 一 个 拓 转 癌 题 。 但 这 种 思想 当时 并 未 引起 重视 。 在 本 上 关 这 一 
论文 发 表 后 大 约 十 年 ， 才 有 人 做 类 侯 的 工作 。 许 多 人 《包括 工程 
师 和 物理 学 家 》 在 不 加 的 领域 中 独立 地 发 展 了 有 限 元 法 宙 。 到 了 
60 年 代 ， 有 有限 元 法 正式 登台 名 。 时 至 今日 ， 它 的 应 用 范 团 ER 
广泛 ， 理 论 上 也 还 在 发 展 。 对 今天 的 工程 技术 人 员 来 说 。“ 有 限 
元 法 ” 己 成 为 必 备 的 工具 了 。 在 本 章 中 我 们 先 扼 要 地 谈 一 下 旦 区 
法 《从 数学 和 角度 看 里 夏 法 与 有 眼 死 法 的 指导 愿 想 极为 丰 似 ) ， 然 
JHdr £88 Rik. 


G, 查 有 限 元 法 创始 人 的 名 单 中 ,还 寇 当 提 到 工程 师 Arg7risiTuatrnerykiartian 
及 Topp 等 人 人 人， 还 可 所 到 物理 学 家 Pra gtr Syuge E McHahon, 这 就 像 一 种 
"onfpypemmbpMgE-H, d mEUrnieueugSud$ EDAM: 

站” 中 国 数 学 家 温 康 在 1965 症 首次 证 明了 “REJ” kat, ix m REA 
可 
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$1 MHHTE i 


1. 从 有 具体 问题 说 起 
HF TRR A 80 7; fio ELIT 28 


| La -- (Pc S E r(zx)ruCm)m für) 


1 (ac z« b) (39) 
Iu(a)— 0 —uCb) 
Eite HRS POQ)€CC" [a,b], reje CCa, 53 
P(r)ym»0, r(r)m0Hr(ioYy RS, 
我 们 把 “* 求 函数 u(rT}yE c Ea, bl BH WW AECSS) SE ERI C390, 
By i A 
Qtstzy)- | C-Pen ta) u(z) e r(r)ut (c) 
Bü —92f(r)u(r)Wr CAO) 
Mule) Eca D BARACCO R X, 
现在 我 们 来 证 明 : Hu Cey T Cro iE 
SOC ue EEpMSOUu DI. ROI, GEocuxcooa.bdE 
OLur REME Q Ce r3, Me ro RmpmEOo Bf, D 
iE: 10 设 aafz) 是 (39)? 的 解 。 任 取 玉 (zcECirap]， 显 
shy =u tF EECa, dhn H 
QtV x) l'E- «Pr tr ert nf Jdr 
--PVF | +È CPV Aryf V MS 
但 Via}= 0 =b) 
Bir EA 
Qt Gy {CPV ter -afV ld Un 
(BRANE, b GD AARD X 
O 将 所 有 在 Ca. b) ESL EGER GRETO EME NC 的 集合 记 作 


Ciro, b), 
Cite, b3s (fir ft€C?ca, bh fins 0-ftbi) , 
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Ld 5685 -* 
OLF (xj) Qf u,Cr - W'OxY] 
=È CP ety e etus WW» 
—2f(Cu, - Wy]dx 
=È Pugh ruzt—2fu, ldz 
+ 2 È CP ugi’ + ru Wda 
— 9 ( Waz 
+È CPW + rd 
但 是 
(EPusW’ ers Mx 
-Puw| + Puly Wa ruWJdz 
-F't-(O uper Wdz 
=È fwaz 
所 以 
OFF (xz) 


e 'CPutte ruo 2f un 十 {CPW s+ ri jdr 


= Qu z+ CPW + rW*1dz 
BH 


Qtu - WJ=QCu) +È CPW tr dr UD 


显然 ， 《42》 RIEAN EC, DIR. E ERWA 
HAF., M 
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- o7 LJ 
h 
| ceo tr) E) Jdr > 9 


实际 上 若 | CPGOW*() ero Wada, RA 
PIW ejt r(xyPF*()RLo,b3EBgdEfaxksksn ME, ufam 
PC) tr) trtr Wr) g 
qdPCIWU(z)—0 
但 P(Gx)20 Æla, blm, HILL 
Ww = p 
pi (rj) 
Bir Ei 
W(x)-C (HRO 
X Wia)c0o —Wi5) WE W)smo 
XU EUG, HEV (rR ET ueh H 
QC' (3) - QU GO) +) PW? + a x QCa v G3 


H R.zE BE QU uC) 13EC39) B ftu oO E35 Elk] f ura 
2) RAQC] u(x)€ Ciba,b3 上 达到 最 小 值 
Qiu,Cr)., WE, [fEHUFCi)€ C$(a,b3, HERRAD, 


BRE 
Qu,Cr) 3 4H (r)) 
=Q luts) PW rW ) dz 
F 2 CP OW) ru (QW) 
-JCW 9d zz Qux) 
Huir) AH CEE, 


[EPOW + AWe +P chan 
tru AW) fi )3dx 


© dxHEfIucD E Cila bit, S OLu Go zz Qf us GD 


H 4 46-12. http://vib.hit.edu.cn/vibbbs 
a 5683 * 
可 视 为 4 的 函数 ， 和 将 它 记 为 94)， 于 是 
PAEAS PW rd 


二 25 CPW eru — JWIdz 
z= Ü 
XPPILGIÉSURAÁALM-BAITCaEgE, XE 
Aat BA CID 
NHEBRISSÉEEE SL, MWE ACO H, ZSC-—o0Np Bo, 
FTES = O0. iH dE nj Ub 
J Gu enu ~ HP )d o 
所 以 
[| Co Ow! erus o Ff de oo 
XE IH EC La, oR 3E. HHEH ELEA 
— (P aia)" 十 rüy—Jj-—0 

qbus x) HR C322 89A, 

S3, AE Ela Oh A 

GE TIB LE) 一 下 全 0 

HMA AAN CEREK TE, EH R dx XD 
o, Hu(x)eC$la,5], dox 

"h 

| Pz) reuke fined 
RAS EF., AË, A 

— tro — FLE) ug ro dEr 

wE E. up. 


" SRAL S poui dz Pup |t- f Puno’ Waz 
= -f ant td 


D &gkdmo—uiizgi£zmbE «Nuus. 
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—GOPCz)usCr)) -r(Czu,(x) — f(r)20 

这 就 是 说 ，amz) 是 (397 的 解 。 

我 们 证 明了 ;wolzx) 是 (39) 的 解 的 充分 且 必 要 的 条 件 Ha 
PROC uCx) Eu, Cc) EREE. 

REH: cRuQ()€C$La.023, 使 泛 函 440) 取 晤 小 值 与 解 
问题 (39) 是 等 价 的 。 

由 于 这 个 发 现 ,事情 发 生 了 逆转 ,刚才 提出 的 任务 是 解 (39)， 
iHe iXokcRu.C x), WERAQ Iu. r)0ER EP. 

里 效法 是 解决 这 个 问题 的 一 御 近 似 方 法 ， 大 致 如 下 ， 

AXoRsQ)eC$Ce. b), EZA 


QC«cz)1- [C CP Gu! G^ iG) rezu (r) 
—2f(z)u(x)]dz 
REMi, dRÁIDGYEEGEBEDAE 
pir), Q,(r), PUT) Q(r)cC;[a,b53j 
( f z:1,2,*"7, n) 
设 它们 是 组 性 无 关 的 。 令 


oae G)muG)0 


i=l 


Gaa a ERRATE ACn OR Tajā, ,G0. 
Has XE 
为 简化 书号 方式 ， 记 积分 


[feoscodi- 0.0, 
把 它 称 为 了 和 9 这 两 个 函数 的 内 积 。 于 是 
Q(u(x)] =f c- (Pu')'urrut-2ful)dr 


(D RAEM, Ea Pair) tapin Rekal.Onm0Xrcca.boüigRir, 
MaRa = 0 (i-1,2,7,8) 
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—OLut.u)—2(f,.u) 


由 上 #4 的 定义 可 知 
Libin +t bu b, Lu, c b Lu DHEA O, OW s A 
此 ， 我 们 说 算 子 工 是 线性 的 。 X 
(f b git bag =b ngt bf 
(b fit bafa, g= bd, 0) + Olf 9) 
所 以 
Q[u,(r)]—(Lu, ui —2Cf.u,) 


E s p. Sap 


Pem 1] 


一 0 


一 S a.a;( L9, e, DRR 


F=, 


E PE 
Eilt, FCa,.,a ER QCu.(r) AE a,, a BERE BEI i n, 
MERR a a, EQ S] 2 FCa,. a, 0) EAM. 
RACE 2r pr BS JIHECKIRIB. Eph 


åF 
a.T’ (i =],2,=,%) 


得 SaLe, p =X, p) 
i=l 


q Lug z (322 x. 


Q 注意 Y ac; 9, 9e0sc To. 0 sacas LP Poe 
LE RE 
*a2,a,CL d mus eg.a CLP uPo Y aia tL. P) 
9 gig CL. re a.a CLP p 
taati i 
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+ YI e 


所 以 Fal(Lo, Poc, P) Ci —1,':7,5) 


(43) 
这 是 关于 ab ea 的 线性 方程 组 。 解 出 (43)， 得 到 au 


42 ;4, 的 值 ， 也 就 得 到 4.(5)-- 入 4.9,(z)。 我 们 利用 
(e, po =) Pr) D+ ry P dz 
--Peto;| +È (Pop, ero m Jda 
-'CPetet irop, Jda 
GEERU,(a)— 0 —9,C5)) 


Biin 
[cover prep jire Bg, ,9,) — h, 
Xx 
Cf ,, ) —b, 
于 是 方程 组 (43) 可 表 为 


可 
1 
1 


| Raat k a+ Ra =b, 


ln A 
| A,.0, 4 LITE T. 4- Rada = 5, 


(445 
l} ha, Rods 十 和 R, 0,75 b, 
或 用 矩阵 记号 ， 写 成 
K-.d-H 
在 此 
K —(RÀ,,] 
47 —[2,,0,,7,a,) B? —[5,,5,,--,5,]p 


© ATRAER ARa MECTHER, Tarts., 
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PVJpupHIABERA ENRETE TEJE ie 《初学 
& ub ARREA TILE P. Xd ei, D SkA 
HER) o 
因为 
R=) CPotot roo ds hy 
BrELSBPERA EARRA. HEHHE REK acen, HD 
人 Ce 


时 ， 
CTKC-0 (注意 CC? 到 CC 是 一 个 实数 》 
(el a 
CKC =p ere Ck | i I= No Nue 
| "ES " 
= > pb; RUE ` c, c, BP, Pj) 
= > e.c, | CPPP, e ro Jda 
5 m r " f 
-PO Xe) (Ne, 
十 (Seo (X eende 
-as( Ye， Eee) 
= B Cur), uT) 
在 此 处 


D EIEEE (0,0,0) 。 
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m" 


— a " 
" Cr)-. N ' pEr} 
—À 


Meg CHIC, e; nua cum HB 这 时 
4.Cx) MB SET. BpeL RELIES D> 0 fiU. di 


c'kc-| Pini c r( udea] nida 
Iri n 


fl dt 
f uds 


可 4n 
CRC sm 
eB FRIES. 
WAIE REE S METE. HO 
| Cd nl) p EC a BS c a 


X 
En 


Mom 


[owo "m zy F'ocoa. V [o z)dz (48) 


所 以 


| d= | {p und de e 


| 
fron ny yas 


eoo) (re 


(heri . 
AF aAa 


ai EXE Aa xL, Parna EEF, 
号 auisitgpo5 S RODIE RA e TR ASAR ER iA REID 4 Tn Mei 


rode o^ A Du = A -o cd o. quz 
Ey paei Da oopa, 
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€ Td + 
Td 
h n : 
[adeno oy woods 
5 
| [Pu'2-ru£]dr 
b 
>f C min Fixzjen trul dr 
© zE En hJ 
b : 2 ， 
2j t, min, PGo cgi ritas 
h i 2 
>f (nincs) e Jida 
=f uf(xodz 
在 此 
i 2 
= am It Lo 
m e h—tu 


HEREKE, HAF RAUDA E E a0, 025, 


Bj XE FG, 0,) 二 M u djR,, -2N aub Xa, a, 
15 =] Em} 
求 导 两 次 ， 得 
oa 
Dear -2R,, 
DE i 
SUP 
aod, ] 
WREE. BIELBEEE P Qa, 0, MEET, Lat ERE 极 小 值 。 
但 六 Car,9,) 是 二 次 图 数 ， 地 以 此 极 小 值 也 就 是 玉 的 最 小 值 。 
MER, Hon NOE GS, De QU uM B BP AZ 
&OtuCG 4 utx)eCi[a.b BIB a. Ma, TETÉSR EHE 
T-TE jt, WEE YTT”, GEJ18 LxhBxpik EZ 
法 )》 求 出 (af eai) 得 到 
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u(z)— X ate, Gr 


iE Edi OE aA” 已 ， 把 这 种 

uir) uar), u E) 
EIBAR AE A REESI” , e Ea pm y s 
Eul) 4 

f |u(x)-—u.(mModr—0 CG noo) 


ti. HE gp EIE DEG DOEA F d SUE 
pun i VLC pop DNCOM 
Pe EE AME OU 3-3UE ESTEE 28 REOR TM 
SERIE CHEE T £4. BARITO MH YEAR 
EPa UDES, H EE ROS L8 $t. A ri 
和 紧张 气质 ， 还 是 来 在- 个 例子 到 。 
B cResÜEEuCe)eC'Qo,1, EERE 
wh a 工艺 (0,1) 
ut0)= 0 —u(1) 
RPR EGL PAP ie) dpmr(r) f(ix)esc-—zgdü 得。 此 方 
TiBS P JE VLRCAR Eu COO = 0 — v Ci BU FÉ V SEE 
u(r)=-r (67) rO 
HEREA RHA, UES, 
HO.Cr)-—2(1—20), P= 0-2), Qua yegtrie91) 


显然 $e.€00— 0 =F,(1) 
He,CcofELa,5 TH LWpESESRE qnd i 
p.(cye Ci[a,5] (i1 —1,2,3) 


D mem “Euin e Cile, di, EEAO am Aea" E a "oup 4 dm] 
i" r1 
Q ut uBLagrange E FE BGGOKkid.oBimug5GtSs. 
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还 本 以 验 证 Pie P ERETTA, MAEA 可 以 选 为 基 范 
数 。 


用 
t 
bum | CPP GO e 9 G9, eda 
(5,7 —1,2,3) 
可 得 
&,,—0.3687 Àb.,- 0.1833 a= —0.1763 
&, — A, B,,—0.1333 Ris = —0.8940 
A. 一 kis A = Bos Rs — 0.0897 C Hi sr UL RO 
X. 
1 L 1 
b=f jp.dzr=| —zi(1—r)dr--—— 
J 0 12 
b- f fpd f s d 1 
"n n 
b f dral amiet d l 
一】 T, r=] -x0-2) + 二 一 30 
由 此 得 到 线性 方程 组 
| 0.36672, 0.18332, — 0.1763a,— — 
| 0.18334, 0.19332, —0. 8940, = =; 
$ 
| —0.1763a, —0. 89400; + 0.08972, = — — 
组 此 方程 组 得 | * 
a = —0.4425, 3,—0.4638, 2,—0.0343 
由 此 得 到 近似 解 . 


u(r) = —0.4425z(1— r)-0.4635c7(1— 2) 
T0.0343r'(1— x) 


T ACD, CLTC PRO, M -Carte iC- Cart Cx 
*C,rzmüTi-C,-0-2C,-C,-C,.-C,-C,, MiC- = 0-Ca.sC, 
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一 —0.0843r*—0.42922z* -0.90602* —0. 44257 


Mubc-l. MRAR OE ul). xS 
uo) EER DORGERDUEHIAXSL AZ. 


jc 2 g^734..-l 
iG eI ee 0-3 


0.90601.) - o. 442s( >) 
— —0.0505 
1 2 3 4 5 6 7 
MERTA =g B" a? S" a? a7 a? 分 别 求 出 近 


代位 wat x) REB Eu CIT )03P 9I E AUT: 


4*1 
"n 点 HT a [E uir) i AR uie 


— 0,0181 -0,0420 


~d 035) 


—0.060B 
- 0,0481 -D.0618 


— ü, 0566 -0.0595 


—ġ. 0504 


-0.0317 -0.0017 


TPXARGBXRAED]s HEEE, GX4DEUEO E AERE =3) 
mjt 
O 这 里 的 准确 值 当然 有 会 人 识 革 
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TEC 

2. 时 兹 法 的 计算 步 又 

上 面 已 经 许 细 地 讲 了 一 个 特 岗 ， 现 在 把 里 族 法 的 计算 步 嗓 总 
2n rr. 

MESTEREK., GRECE OM rn m EU Edu (0039 满 是 给 
4E ES XL ELA PER EBE Ci xp. MoCila6d. $E 
Phoe, MupHRHEdpBDUSGD o RI n.0r2c€M, W iG 
QiuCrc) BU Bb, REA AE 

TUB E EO db RB Ar [RI REL ZE TRE, P Je 

12 (E IR EE JS M 中 . 选取 线性 无 关 的 几 个 基 转 数 Pi ,ps 
Pae IF 


uH. x) a Eyy 
ml 
向 


沁 Msi Sapo ee 下 一切 实数] 


1 =, 


IAI iZ 
Qu E= (Dav tr} ) 


ef d.i :9,) 
GR AG. C RGBjgDe.4..0. Hn JUR EE, 
2) sRHjat,2a2,-,012, WARF (Cai. o TE Cat, a9, e 
as) 处 取 节 小 值 ， 即 
FCat,ob,-.gi)-— min  Qiutr) XD 
di UR OL = o 


1335J|£etEZ; EHA), BJ 


,d= 
解 此 方程 组 ， 求 出 (aĵ, a}, e A) * 在 一 定 条 柱 下 ， u,(r)- 


(QD B4 min Qtun ma Quiz] 
“T EN, = Tii 
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S ato, ir EEA u GEM , rp RM 


L=] 


QUu,(r)j-— um. Qrutcr)), B 


f Iur) u 0 M n) 


(uolX) 是 变 分 河 题 的 精 闯 解 ) 
这 融 是 星 兹 法 的 基本 轮 记 。 


$2 AREH 


fk 9 E BOR US ES OLu Ge TRO ERICH CD n C ) FRE ADLER 
f fic XE Ec SEE SEP (T), ,9. Cx) C(BBLESCEEPRLBI SESS 数 )， 


Su r)-— Yapa), EAC ONER A E 28 ML E 


r=] 


HORRANB. EXE u CEA CORR. RAA dE 的 指导 
忠 想 是 容易 理解 的 。 一 般 说 来 ， 也 收 钱 得 比较 快 《 指 n 不 太 大 时 ， 
吗 莽 比较 小 ) 。 但 是 ， 里 兹 法 对 坛 探 函 数 p,(z) 的 光滑 性 要 求 较 
高。 对 于 边界 曲线 比较 复杂 的 二 维 区 域 (高 于 二 维 的 区 域 碗 不 符 
二 ) ， 蓝 梅 造 出 这 种 试探 性 函数 往往 是 困难 的 ， 有 时 还 是 不 可 能 
的 。 玫 外 当 半 半天 时 ， 夭 度 提高 很 慢 〈 注 意 这 与 上 文中 所 说 的 
"36 n 不 太 大 时 误差 比较 小 ”并 不 牙 秆 》。 面 对 这 种 情况 ， 有 些 
人 对 PAT) 的 光 油 性 隆 惰 要 求 ， 干 脆 选 取 分 片 多 MA WFO 
作为 9,(z)。 这 种 基 畏 数 仅 仅 连续 ， 在 有 些 点 连 导 数 也 没有， 这 
就 扩大 了 试探 函数 类 的 范围 。 当 时 ， 人 大 们 难以 佑 计 到 的 是 :这 种 
做 法 居然 引起 了 计算 方法 前 巨 天 变 草 ! 产生 了 现代 的 适用 于 广阔 
仿 域 的 有 限 元 法 。 本 节 打 算 从 变 分 法 的 角度 介绍 一 下 有 有限 元 法 的 
指导 思想 。 

再 来 看 一 下 问题 《39》。 我 们 已 经 若 道 ， 只 要 求 出 uw), 


T BEER uin gun EE HR, 
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使 (40) 式 所 岩 示 的 泛 函 QEu(z)] 取 最 小 值 ， 则 ule) E (39) 
的 解 。 但 泛 函 QEu(x)] 可 接 下 式 进行 变 搞 ， 
Qute) J= | C- C^ G)u! G0 uz) 


Tru'(x)—-2f(r)u(x) dr 


——PFí(x)u'(r)u( x) 
«ccs Lrut(r)—2f()u(2))0dx 


-l'tPurta ru—-2fuj]dxr 


《注意 #(2)=0= nu(b)) 
现在 ， 我 们 用 有 限 元 法 来 求 xe(z)， 使 泛 函 


QC G2 3 UP u^ tru - 2f udz (46) 


Eul EH EARN (B 
1) 选取 分 片 多 项 式 作 为 基 郴 数 〈 这 是 关键 问题 ) 。 我 们 要 
JEB —^7 H8 22 Er RETI LE BC. Ca), EG TECO, b JERE 
Vix )-—u(z,)-—u, G—0,1,2, 7, 
(tar, bz 
CD, Eg 13? 
Ry zr 
取 Vu(r)m 一 — 
Lo H r ELE 2] 
(r-—zx)/h, Mx Etro x) 
Q(r)j— (Ts r) h, HIELT, Te] 
Co H or Elta za 


二 ”为 便于 推导 有 关公 式 , 我 们 在 原来 的 Qu GO 3 表达 式 前 第 内 地 。 Ziu.tr) dE 
(itPutterut nf oodz GR uM Iu co bit f CPU tarut- 
2fuldz Rahi. BDOIIGE. LE SD EU 
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& 81 æ 
(£= Xin R; MW rC€[r,..m,) 
Qr)- PIN C 
0 Hrer, JEE] 
X 
WT) = j 0 MrE€[r,.,] 


{r= r, h SIENTE. oa. 
Pt} 分别 如 图 135、c、g 所 示 。 


ES zz A2 qu 
1 | 
] J 
上 1 j 
| 1 

l Fe (a) 

A. RI I E 
| d l 
| 1 | 
| 1 Í 
|! — mn | 
VS 
| 

aeee a 

y FT i 
| 
! l 
| 
! 

l Í 
1 | 
9. xA 
| 2 
E 2 E 

at A 
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= B2. 
TYrli3e«. b. c. gd 联系 想来 观察 ， 可 以 看 由 


V(r)= S n,o, r) 
SIZFYEERITSO Rr), €.(r), e PRETAS NES ERG UR 
情形 〈 即 分 片 钱 性 和 多项式》， 在 有 限 元 法 中 ， 这 是 很 重要 的 基 郴 
En 现在 将 有 限 元 法 文献 中 关于 %,(z) 等 的 通用 术语 介绍 一 下 : 
区 [El NEN 称 “ 单 元 ? odd fE Eio xà X "3 zr, 都 
m] €". NEXTAR 
NED Cxr,—r)Mh, 


(47) 
nr, 
PEME "TERT, 于是， 下 述 基 质数 可 表 为 : 
(QN PP Gn) Us y Cg, 
"NOSE ] ] n (48) 
XU oo«j«n. 有 
[VM I M rCe, 
py) NDT) ard 4 Eu (49) 
| D uw r€e;je,, 
x 
Ngn {r} yr če, 
e r— S (505 


‘oo Mr Ee, 
这 些 趟 子 虽 然 难 展 些 ， 但 只 要 用 得 史 了 ， WRAT. 
RE u,—P,, HEt, Uap ttt. u [EDO D. mr eg dr II 6r) 


= > ug.) HAEIRIERHNURV. 用 集合 记号 表示 ， 就 


下 


是 
Me FS unco Duo Pons ae enc T 
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IER DEM,’ ERU (a) P, 又 内 为 
0 Hi -:D 


y 
pau) 1 

o fre] 
Br 以 


Fia ag Gn), mu, 


XX RE, dRÍUDELÉERG TP 
aSr rae Iro xb 
HEE a, D 划分 为 Fr 个 单元 e =r, aar 《7 一 1 二 )》 
Xx T EEHB BEY T) €.) -, Peh, FEE T XXE X 


Vix)- Mu) 


2) 按 单 元 写 出 单元 矩阵 
利 星 兹 法 一 样 ， 下 面 的 工作 是 要 在 试探 嘎 数 业 上调 ,中 ， 求 出 


LARR uu (n) 人 pz ERL Cr) EX uu Cz) 
i-4 


上 取 最 小 俯 ， 即 
min QF (r]-2 Q[u,(r)3 


VESTEM, 


HVEM, ih 
Qula SQ D nie.) )e FO su) 


注意 到 xu 一己 ,已 经 事先 固定 ， Bs pL 25V (Cz) EM ,类 中 变化 时 ， is 


(pud D Eje. py), 9. 0r) RATE, LI: NEM PO 
i= i] 


=0 , Ære T,, Rj* leget) me, 0 (220,1, *", n) 
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WOC (CE). u. c, 2 的 ?元 函数 。 只 要 再 用 《这 也 和 里 
效法 一 棕 ) 

o 0 €i —1,2,-,9» (515 
FERTH ta u MERETE, HRe ut ut 就 得 
到 近似 解 


un) > ute) Oud s — PO 
BERRE, DIBIENE DAE ARTE, MIK ERRA 
BATEK, MASRAM, ZR SUDLDEEDL CD B Rd 
说 来 ， 因 为 


QU JSE CPV ert 2i Ms 
-» [i CPV rV? 22fV da 
=] "jl 


vli CPV’ tr Vf dr =Q, 
s A rj dd 


-=l "ld 


在 此 ， 
Q,-1| [PV 5 rV:-2/V]dz 
à J*, 
注音 V(r)= Suo, io. T IET 
s eS e, 
Ee; 单元 ， Xj 之 TT 这 时 内 "ob | + S2 
有 PX) H9; 0) RAE OO Aj Aj Ajri 
图 14) 。 再 说 详细 点 ， 福 8; 单 m 
Ju. I 


Q,(xr)-— NP T= th 
9; (r)—NfP"(r)o-(r,—x)h, 
Pr EE dg e,.E 
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sps 
V(x)-u;,9,(r)-4u PET) 
—u, VN IP rtu, NiO C) 


CPV” tryta] 


dN dNÍ£^ 
Lr E Ue) 


priu N £P pu NN £P Y -2 f Cu, IN EP ru IN TP yd 
I ZU G( EE EH "s FON £P » Jaz 


AF CP AFP 


1,2 | ANN ren :| 
i PP (P( d ) +N £P Y ldz 


l 


=i, fN P de= u JN £P? (x) dz 
i £f 


Ap 
u ANI UP d N G C "T "T z 
vue], (rt Pot Je 
(52) 
sief, rode 
T 
(2,8 =1,2 P—i.en) 
ja 
Q; = kju} TE IPUIM My — sisi bius c biu, 
我 们 要 写 出 导数 
oQ | i òQ; 
T 
来 。 因 为 
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—biu, 


i+Í 


二 —biu 
Pr vl 
Qi lkitiubadkiiin.u., 
th —bi*'u, zi bitu 
WUUBHERRO BRE. BL 


9Q s» 30, 9Q., 90... 
Qu, = ðu, | Ou, du 


—hiu.,. tÁ u, b+ hitu + RItiu, bit 


现在 ， 只 要 令 守 一 一 0， 就 得 到 关于 ”的 线性 方程 组 了 。 但 是 ， 
我 们 打算 在 每 个 单元 e, PRERE, AHA nm SA Ru PE. BT 
以 ， 对 于 Q, 而 言 ， 我 们 先 求 GO ASIE 
i -= kit, +j, bi 
(53) 
Pt hb 
EHEHE. $ 
| [5i 
i= È | = B! 
[kii ki; lb; 
于 是 
aQ, | 
Qu - T 
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K'XBSOUBGLe.HWy “RICER” o Rmus 
元 6 VEBABU ER. RAH ETETE ERK UBLDU' 来 ， 就 完成 
了 第 二 个 步骤 。 

3》 和 将 刚才 得 到 的 单元 矩阵 人 ' 和 2#'， 组 合成 总 矩阵 ， 写 出 
线性 方程 组 。 

这 一 步 又 在 许多 蔬 上 都 是 通过 例题 来 说 明 组 合 规则 。 但 只 死 
idiH e BLUR. Pype BGUIARBS. TERRAE HE 


HEXIA, RITESH S HRE 
ag | - " 
du. = 0 Qi -1,2,:c,.,) 
也 就 是 
OQ AQ r 
àu, Qu, T?) 
9Q. ac , 
ðu, "Ou, - | 
0 Q, Qa í 
Ou, Ou 7. 
E | 
| 
a Q, N 
Ou, TUE 
用 2 A, 98 GEBRi,-RID 
Sd Aue NETT. 
ao. 
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ðQ: 


— h2? 2 Z 
= kist, + hit 5j 


Óu, 

Sei ki s e khub 

de boue idus — 6 

Pics Hi dath haha 
e =R] H, at Aiz Hp 0I 
一 Ai HU, Sd RI M. .—Ó5.! 
5 =kiu, Riu, 51 

e. = RI a Ral.— bI 


将 它们 代入 前 商 的 方程 组 ， 得 


(ith? u tA? utelu bl—b!=0 | 
i H ALI tÁAlhu,—-5$—5j-0 


ter 


ko Haa HRR PRT Enas hi;'u, | 


Kc 


—byi.pbt'IÍg ) (54) 
kelu, a RIP AAT u, hl, | 

-bpi-br-0 | 
Riau, (d R$,u,— 05 —0 j 


EE JA AG RJE 


lbbbs 


http://vib.hit.edy.c 


H h 16-14 


pi Mi nr 


=M 


TAS ARRNE ERE RR RAORNS RRÀ noB RREORRR RAE ERE NARRA RU da 
1 
1 


性 d T*iy ir 
0 Cip yta o i 
0 二 
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bat bi— kitin | 
| biebi | 
| bibi | 
P= : (56) 
| b;trbyt 
| bii rb | 
b 
U?-[u, ,u,01J (57) 
WARRE 
KU = 8 
我 们 设 u(r,)—u,73 EAR, Bik PRE PERPE BOE: e E 
组 (572. 


in€urQj =u Aulroo osun R ERM, g5kAmt s n 
uuu, 了。 而 是 在 列 方程 组 时。 C. = 0 这 一 方程 也 失去 
AXT., TE, fEXXRRORPEBRORÍIGIBEUIS TOEAnhEu, uy 6, 
2 1 的 方程 组 

(kitki s +e u, = b+ bi kiat | 
kiu tiki Ri uw, tki u= bib] | 
kiju, s CRI RI Duna t EIU. | (58) 
sbl po] 
Raiu, a ERII + RI) He 
=bg bik, 
HRO(58) x66 AB BTE 3. 
K*U =B (59) 


在 此 
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f bi bli—ki., 
bi d bi 
bibi -B (60) 

bijrbl 

b Hbk. 


又 总 刚度 矩阵 为 
(Rist+ ki, Rh Dag 
| hii? hip B hls RIP 
i "TE kirl + Ai. 
(61) 


观察 (505) R GD 式 ， 可 得 总 刚度 目 阵 乒 和 兵 * 的 “ 拼 污 法 ， 
这 是 进行 腿 元 法 实际 计算 寺 必 户 知 道 的 方 法 。 现 在 ， 对 矩阵 
六 * 的 “ 捞 谈 法 ”说 明 如 下 GERK nu DEBEAS IK ^89 Ut ae TR TH 
位， 仅仅 在 矩阵 的 右 王 前 有 些 不 同 。 请 读者 自行 分 析 ， 最 好 结合 
算 例 ， 不 难 作 出 有 关 结 论 )， 

对 边界 单元 e;， 算 出 i。 对 单元 24, 算 出 算 阵 


ki. kia | z 
^ "D k? e GRK "为 单元 ex 的 贡献 矩阵 ) 
21 zt 


| bes 


—R x, ADAGecEW, HEARR, 


hl. Ris 


K= " (j—2,8,7*,5—1) 
Rta kis | 


(注意 Ri Ski ARRE, ABD, 

再 将 这 些 和 矩阵 乒 : 按 对 角 元 素 位 置 “ 焊 接 ” 起 来 ， 就 得 到 抵 
阵 K*。 焊 接 方法 如 下 表 所 示 ， 
AI, À 


T EFRR K, — 


i rn Bi j 


//vib.hit.edu.cn/vibbbs 


htt 


» 02 + 


H h 416-14 


see. 和 
PET zT E: 
t u E] + 1> u T RET 


DLE 


f = 加 


agikem'tU qug uH" TANTA 


-" 
号 


rH 
H 
H 
daRrhraa T EYETWALARAATASRERARE RARI Lau RR RA 


| TRE Aeulaaara BER RlldARq4ARRTRGSQIRRRRR Gi ugqus punk ddd db up kr RH URRRS RAE GRERRERRUN ea k d E Eaa ehh E aaa 


--- - ASR-TARRARLATAWARRIRSS TRAY Dul-4un-34 AT AR dou 


' M : i i 
i oy +R Ki : 
l : -aRh"RRGRERad-AGiyrAmTRR-LETARTRuan-nuuua bm mmu LLLLI 
| i ; : 
E in t! E T E 
19) ， i :" ] g 
d | : "iy ia tos i EI : 
i ---:1-9 Xl Ld cnl ECC CC DL ODE sub CM CR IG t 
E TI HI T i T. : 
: EE ie} tot i :H i 


: f f rt : L 

E i +y i i 

: | Bes TUM uo: 

DUE EN Pr d Fhaddr H 

| : 

: P 11 tt : 

i H E Stud 
~ 3 : 
| ————————  Pw— do 
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. 03. 


金 于 矩阵 如 的 组 合 ， 可 先 接 单元 8 / 求 出 矩阵 


Bi=| ’! 
(n 


BRE C60) XE B. 

4) 解 线性 方程 组 《57) xk (59 。 也 子 电 子 计 算 机 正在 普 
上 用， 大和 何 解 线性 方程 组 在 这 里 就 不 必 多 说 了 。 只 说 明 总 加 度 和 拭 阵 
际 了 具有 对 称 性 、 正 定性 之 外 ， 还 具有 区 朴 性 。 征 和 阵 的 每 一 行 至 
多 有 三 个 非 零 元 素 忌 ， 这 就 大 大 践 少 了 计算 量 。 这 一 点 也 显示 了 
有 限 元 算法 的 优越 性 ， 

用 有 限 元 法 能 问题 《39) MHA PRCE, MEY 
一 例 ， 详 细 说 明 计 算 过 程 。 只 要 把 这 个 例题 的 计算 过 程 全 部 弄 异 
了 ， 也 就 进入 有 限 元 法 的 太 门 了 。 

例 求解 


和 rc 0, 1) 
u(0) 02 u(1) 


《注意 这 就 是 在 上 节 中 用 里 兹 法 解 过 的 例题 ) 

我 何在 讲 有 限 元 法 时 ， 从 始 至 终 都 以 变 分 法 为 指导 和 起 想 《 实 
际 上 有 限 元 法 正 是 在 变 分 法 这 棵 二 树 上 开放 的 草花 》。 为 了 使 初 
学 者 对 有 限 元 法 的 指导 思想 更 加 了 了解， 我 们 在 解 这 个 例题 时 ， 先 
不 按 上 文中 兽 绪 出 的 四 个 步骤 来 计算 ， 而 是 更 直 返 地 从 泛 量 写 出 
线性 方程 组 ， 再 求 其 解 。 

我 们 已 经 知道 ， 只 蛮 求 出 wkz)E cilo], BE TE ER] Qu C )) 


=| Cuca) Tut(r)42ru(zx)]dxERdR MB, Bu.Cr)dkEPUES 


的 解 。 用 节点 0 H qo Š, LAE MICO ,1 J 划分 为 4 个 单元 ， 
Bp: 
D ”起 刚度 第 阵 的 每 一 行 欧 堆 元 分 布 ， 汪 单元 的 身 号 ， 排 列 有 关 。 这 一 点 ， 在 多 
FA EIJER E. He MEA RAEE. 
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1 3 
CEEE 1) on ( 3) 
£. -G : | ) 


az 一 D u, p(x) 


= 


Bu = 0 su BARE, BLU 


4 


u,(x)— S up x) 


在 单元 e,= ( 0。.-， 及 其 两 端点 ， 显 然 试探 函数 u(r) 是 满足 
ui( 0)= 0 及 w,( 汪 = 的 线性 疯 数 。 由 图 15 可 知 ， 当 zx EL, 
1 

LP 


ri 
ux) dur 


4 
ut(r)-—4u, 
ia 
a. Fig 
t, udi E 
l i X n t 
j dio 
"d “a t 
e zs M 
4 2 
图 15 HE 15 


TET- asfi i). 由 图 16 可 出 


nm Manso -1/4 
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* 95 « 
—-u2—4t)dtu,(Ir—-])9 
uj(zx)o4(u,—mu,) 


Mr ee (t JJ. 由 图 17 可 见 


uamue) Etu -)/ t 


4 


= u.(3—4r)- u,(4r—2) 
ui(r)z4(u, — 1i, ) 


Mree,-(5. ! ] 时 ， 由 图 18 可 见 
sr) ul) faz) 
uir) = —4u, 


ER QCuG))- SQ. 


| m j 


Q, = [acuto ul(r)-c-2ru,(r)]dr 
4 


-|. (ui -ruirzruQjdr 


现在 按 单元 算出 ，， 


Q,- sieut Tüi4r)!'ul- 2z(4áx)u,ldzr 


© &3 (45) 式 
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* 06 >+ 
«(4 + iyi " du, 


Q= [tco ann (Az-—1)u,)--16(u,— 9) 
4 
-2r((2—4x)u, + (4x C1), ) Jd x 


1 3 
Mc utu HAU — uY + 35 i a) 


s= [icus Ar) + u,(4r—2)5)!' 16(u,— u} 
-2z((3-— dT)u + (4x —2)u,) dr 


ig t us) HAC ua — u,)* 十 5 (uL uU.) 


C161 — 30* ui c 16u$ b 2x 401 — 2), jde 
1 2 7 
"(t 4 yit 32"? 
But. EUER Q[u,Cx) EXT 2E, Uu, ds 的 二 
USER BR 
29 = 0 (i =1,2,3} 
i 


ELE SOETu, to ts 的 线性 方程 组 。 从 刚才 的 式 子 中 再 次 看 

di Q Eu ESAE, Q Eun s.p XE. Qu, ts 的 图 数 ， 

Au COE Lp EIBO Es, s. AARO, TÆ 
àQ 23Q,,2Q,. 


Ou, du Ou, 一 (5*3) E 


EGO u)- 8Cti,—u,) + 二 


s; 
219 sor (4 - 8 和 + 


他” 为 秒 一 起 见 , 在 计算 定 积分 时 ,一 站 用 馈 形 公 趟 或 币 形 公式 求 近似 值 ， 不 管 能 
TERHREY. MEADES pidr (二 jb- a 


TP 
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+ 97 æ 
3D. ðQ, 00s. Tu sus - 
Qu, Ou, ur zii FH) POM) 
tSt H -Fu4,)—8(u, =u) tz; 
"I7 7 8 
"i e) i 16 二 5 二 8 Jutan 
St ac l (ut ua) 8Cug — 14) 
5 m 
ur ttem 
(rye 
aygo de 一 0 (i 一 1,2,3) 得 线性 方程 组 如 下 : 
16.25u, — 7. 875u, =— 0.125 
=—7. 8708, +16.298,— Y .675H, = —0.28 
—7.875H,-F 16.258; = —0.375 
再 表 为 矩阵 形式 ， 即 
16.25 —7.875 0 | u | —0.125 
—7.875 10.25 一 7.875 |j u, 一 | 一 0.25 
0 -—7.875 16. 25 || us | — 0.375 
也 就 是 
KIU =B 
K*By ^3ARBJEESHER? ， 正 站 上 文 所 指出 的 ， 它 是 一 个 对 WR. iE 
EH Fi pt A EEE 
解 此 线性 方程 组 ， 得 
ut 0.0354 
u$ =— 0,0571 
ut = — 0,0439 


于 是 ， 我 们 用 有 限 元 法 得 到 本 例 中 微分 方程 边 值 问题 的 近似 解 
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= OB s 


3 
ulr) = N utp {ry) 


LS! 


详细 写 出 s.(x) 的 表达 式 ， 就 是 
—0.1416 x 


see(o.1] 
—0.0354(2 —4x) —0.0571(AT —1) 
ase(1. 1) 


u,(r)- — 0,05713 —4:) —0,0439(4r —2) 
1 3 
一 < [ ] 
Mx £ d 


—9.1758(1— x) 
3 ) 
8r e( 5. l 


db. Persa o 0 Loup 


E (i)--* 0177 
u, (1)- -0 0354 
Pp (3 )= -0.0483 
cun (1)- —0.0571 
u, ($)» -o 0505 
n, (3)- -0.0439 
", ($)--o 0220 


TEILS RTE GEEA NARRER CE 20 ， 以 便 观 
察 误差 分 布 状况 (参看 本 章 8 1 中 的 表 13 
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ge 2 


| 节 5 1 4 BEBO CLUB | KREEME 
z= i -0.0191 do lem (|  -emm - 0.0420 -0.0177 
z= —i— -0.0350 — | - 0.0569 — 0.0354 
r= — 0.0484 | — 0,0618 -0.0463 
z=- = 0.0568 | - 0.0505 - 0,0571 
zeo - 0.0578 -0.0327 | - 0.0505 
z= Ê - 0.0508 -0.0142 -0,0439 

r= - 0.0817 | - 0.0017 - 0.0220 


ium ERE. GRRII XXTEHUUSIBIENSQ, BONSORS 
MARETE. MARTER S AER AnA, H 拼 
成 总 刚度 矩阵 玉 (或 KR*) 及 8B 的 方法 。 为 使 读者 能 掌握 这 种 方 
法 【就 是 现在 通用 的 有 限 元 算法 》 ， 我 们 对 于 刚才 这 个 例题 ， 用 
上 文 指出 的 四 个 步 又 理 计 算 一 次 。 


第 一 步 ， Quo. i zi Ss 1, 将 区 间 Co, 13 划 


j ] D. 1 t 3 3 一 
safo 3) (的 G De erna 

第 二 步 ， "n TERETA AER, WAEA E, dU 
按 单 元 计算 时 要 用 的 公式 汇集 如 下 ; 


dNIP d Nj? 
ZINC qi gb + NIP Nin js 


bi=| [NP dz 
"i 


T ERRAT, TREE bs. 
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- 160 = 


NiPDG)-G,;—z)/h, 


Ch —r;-— rj) 


Nn (r)-(r—r, )/h; 
在 单元 2; 上 ， 要 算出 


[Ri ht i i 
Li LN TN b | 


“村 | 
Ri, Ri, E i bi | 


因为 u(t0)=0==w(1)， 所 以 我 们 要 用 C60) 和 (61) pru K*4 
B, 


现在 xz,=j/4， h Cj —1,2,8) 
所 以 


NEO =h- x)/ 1- jac. (IN Er (zy = 
N 1? Gy) (a -4 )/ aury 


EN £D (aI 54 


x 
P(r)m1ur(z) flzj=-r 
所 以 
L 
khim LONED E+N E Jda 
Ta 
l 
4 
=| C26 G - 4x) dz 
UA 
wyte G-no) 000 2a 
Cui 18 
了 
kimki Ef UNAN EYAN P N f 3dz 
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a IOI 
! 
=ff c-eeg-4nüz- Pes 
MEL 
1 63 
dmi T EET. 
i 1 
pls=|, necari 十 1)Jq7ee 415 
4 
x 
L 
bia [^ Cong -andz 
1 | 1-8j 
mr) DE Ar) m ipil m 64 
àvAÀA E! 
1 1—3j 
d 
Pr E 
hl 
ak -5| (u-z es 
K'- (0 Bis 
63 1 | , 
-18 *16| (a 25/)/64 
-$ 
B (60) ERTAN 
1 1 63 
4i6 t tie ~ie 0 
63 1 1 63 
K*= -16 dis * 416 15 
63 1 1 


102 
| 1 63 3 
Be  —18 “i 
~- 88 gl 63 
16 8 16 
69^ uet 
16 & | 
Hd (60) 式 得 
^" ` 
(arroi TOT. 
| 64. 
| 
xp diha 4 19. 
B pi Z | 
| 12 | 
3 4 5. i 
DILE 64 | 
TE, SGDXCPRAUERS, uu, zs 的 线性 方程 组 
l 68 站 
Se 1 dns 
63 1 63 u 2 
-ie êg% 716 “` —38 
63 1 3 
| 0 一 证 zu lis] p 


此 方程 组 显然 与 《61) 相同 。 

第 四 步 ， 解 这 个 方程 组 。 

这 些 就 是 用 有 限 元 法 解 1 维 问题 的 实际 步 又 。 补 学 者 如 把 这 
一 节 的 内 容 掌 握 了 ， 就 不 但 学 会 了 关于 1 维 问 题 有 限 元 法 的 实际 
计算 (特别 是 关于 总 刚度 矩阵 的 拼 普 法 》， 而 且 懂 得 了 为 什么 要 
这 样 算 。 和 否则， 如 果 只 知道 机 械 地 接 规 则 计算 ， 思 到 复杂 的 实际 
问题 就 难以 处 理 了 。 

本 书 是 在 变 分 法 的 共 础 上 来 介绍 有 限 元 法 的 。 试 图 使 读者 在 
懂得 变 分 法 以 后 ， 进 一 步 领会 有 限 元 法 的 数学 思想 。 在 这 种 数学 
思想 指导 下 ， 透 惹 地 掌握 1 维 问题 的 计算 步 酸 。 关 于 多 维 问题 的 


振动 论坛 http://vib.hit.edu.cn/vibbbs 
* 103 * 


具体 计算 及 实际 应 用 ， 在 这 里 就 不 谈 了 。 读 者 在 看 完了 本 IE 
后 ， 应 当 参 阅 其 它 关 于 有 腿 元 法 的 书 ， 使 自己 能 况 在 理解 的 基础 
上 掌握 多 维 《特别 是 2 维 ) 问题 的 计算 力 法 。 

最 后 ， 还 须 指出 ， 现 在 已 经 有 关于 有 限 元 计算 的 多 种 通 用 程 
序 。 读 者 在 熟悉 有 限 元 计算 以 后 ， 再 学 会 使 用 这 些 程序 ， 硕 能生 
电子 计算 机 上 迅速 而 准 勿 地 解决 各 种 实际 癌 题 了 ， 
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第 二 部 分 外 推 法 
第 五 章 外 E. 


外 推 法 是 一 种 简便 易 行 而 精度 又 很 高 的 数值 解法 ， 在 国内 掌 
所 了 这 种 方法 鹿 人 还 不 多 。 但 这 种 方法 实用 价值 很 高 ， 应 涩 介绍 
给 我 国 工 程 技术 界 。 在 本 章 中 对 这 种 算法 及 其 应 用 作 一 较 系统 的 
介绍 。 


$1 外 推 靶 的 实例 


我 们 打算 先 谈 两 个 最 简单 的 问题 ， 从 这 两 个 屋子 中 引出 外 推 
法 的 指导 思想 。 

例 1 求 5 的 近似 值 
我 们 知道 直径 为 1 的 加 的 肉 按 正 # 边 形 的 周 长 为 


„= asint 
我 们 知道 
d, —asinZ f T ea (3 n — c0) 
现在 求 由 dg。 作为 z 的 近似 值 。 


先 取 nn —12, f d,,-125j.—3.1088286 H Hia —24, 


12 
得 d,,—-245in;. — 3,1326285 这 两 个 值 误差 都 比 iz. HAE, GU 
果 我 们 用 下 而 这 个 算式 


d*= da c dy) 
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Xp .bExRd,Kd.,ixbf] *5MEU RE 
d* «(4x 3.1326285 — 3. 1058286) 
—3,.1415618243.1418 
是 不 是 碰巧 这 样 呢 ? TPu—iX. XX m. ERIH 
. Œ . d " 
cda = 6sin = 3, d,,::i2sinp,—3.1058286 
再 用 刚 半 的 办 法 ， 得 
d* 5 3- (4d, d.) —3.1411048 
u,183 8 Had Rd Veni e i x em o C. 
一 般 说 来 ， 我 们 先 求 出 #5 的 近似 值 d, 和 da;.， 再 用 如 下 外 推 公 
式 
d* (4d, —d.) (62) 


i 18 s EE d ud, ES BSEO dt. ix FpHEEDBSSRUE TU f& 
* 外 推 ” 出 较 精 确 近 似 值 的 方法 ， 实 在 术 简 便 了 。 只 要 对 原来 的 
ETEME Md .进行 三 次 算术 运算 就 行 了 ! 这 样 好 的 算 d, 


我 们 在 实际 工作 中 为 什么 不 用 呢 ? 
现在 来 论证 外 推 公式 0625 。 由 熟知 的 Taylor 公 式 得 
sins =a curata ST ue («c c1) 
I 5! 
所 以 
JE 
om " 1 m cos- z5 
SN nT hi 120 ns 
T 1 m? ] 
Spog t ng 
在 此 
eu Cn 
i 120 
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" 106 * 
* 
C Te 
于 是 
mc r- a 1 
.c- HSIN "S ze E EE -L C Lus 
m? 1 1 
人 n? 十 b 
m? 
* 于 二 全 则 
NES 
d, =T + n* (535 


在 〈65) Rh, WORCÓR n DO. WRC HRS 有关， 但 
C. 有 界 ， 即 存在 与 ” 无 关 的 正 数 M = Lo. ft 


所 以 
d = E -— 二 " 7 — 1 : En 
^o 75 (n3 7 (ny 7 74m pent (64) 
于 是 
CES 
4d,,—d,--3z-F 1d e 
1l ET ROT OM 
g (id cds 4 ;( 4 一 Cs 
令 
If Can . 
i( d C. Jess. 
Di 


1 
d'y Gd di) uar (65) 
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由 03 及 Ca» pn 
m 
d, =T — o(*) 


d, O(s) 
RREH, RA ld, 一 zi 和 ds, 一 可 均 不 低 于 去 阶 《 在 许多 文献 


中 说 成 BO zs 阶 ”》 
由 L685》 式 可 知 


d*-—m: ofa) 
这 就是 说 ， dt=Hada o d) RERI PRE, 可 见 d* 
比 d, 和 d;, 的 误差 小 得 多 ! 
$12 求 定 积分 的 近似 值 
对 于 定 积 分 
TUNE 
fUr EERE a a 1 ] 分 为 N 个 
N 
Laura E UE z=, n= ta N= 
siei 


Jis fiu ut Ln eT Tus {Ern fyw 


ER 
LA o + Jupes t =) 
e [foods 


bens MEE ME, GE EH BOE XS SIBI IE(U 


1 
D P. EN m <H (11 与 1 无关} ， wie R= obs) 


振动 论坛 http://vib.hit.edu.cn/vibbbs 
» 108 » 


信 I wn 及 1 x 用 外 推 法 处 理 ， 得 到 更 好 的 近似 值 。 
s= h, XEBA l, HEER GHHA ma pk (Euler 


Maclaurin DAL E AW FR, BD. FARGA yf 
(0, 10 ER2N +20 R (577^ E, Hu] 
I5 AC AACA t CA He t nA NA Cpa AN 
(66) 


1 
Gsis] J(xodz 


CoC Cr CRE àh ARDRE, RACs 有 界 CONI xe 
YES IN JG XA EE M, E IC uu] SM Xp EAN SEO D, 所 
以 
I, 4C A+ OCA) (657) 
(66) EK (67) 式 指 出 1 与 定 积分 精确 值 的 关系 。 它 们 都 称 
为 1 的 新 近 展 开 式 。 对 外 推 法 来 说 ， 关 键 问 题 就 是 求 出 这 种 渐 近 
展开 式 。 有 了 这 种 展开 式 ， 就 可 以 用 外 推 法 求 尘 更 精确 的 结果 。 
做 法 如 下 ， 


BM CYPEE SEU 
， h £ ht 
nyat (4+ oli 
=l ete h*4- OCA*)0 
所 以 


Alas C AHA OLR) 
F 
=47 CU TO?!) 


Luül,-1)-1 +O CHE C67). 式 代 入 此 式 左 端 ) 


这 就 是 说 ， 如 果 用 Dre ul 1o EX I= | /G0dz 


D 关于 Ruler-Maclaurin 公 忒 及 其 证 明 , up BckinOEGDED ADGRDLX IRE 
UI P p.157 160 
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的 近似 值 ， 那 碌 ， 它 的 误 姜 是 不 低 于 如 只 的 。 原 来 用 樟 形 公式 求 
tigre. 和 只 有 不 低 王 户 阶 的 误差 { 见 L67) 式 ) 。 

我 们 已 经 找 出 了 比 梯 形 公式 更 好 的 求 定 积分 近似 信 的 公式 。 
现在 把 求 1* 的 算法 直接 写 出 来 。 因 为 


L= r-o OES COCE (3). (5) 
Aee heran] 
=- [fo 4 (t) 十 (x) 2j (3) 


ea) ser) ro) 


所 以 
1 


I= Tas m |I HG) Air) 
Ge) eO) duc ) 

bd MN Ie 

ad Jes n. m e iO) ] 

p= HM Mr) + ue) 4C) 


CAERE T GOSETE 


-wí fC0) + "era 4 24 AR) 


D Bh of 二 -:0 
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eG) nr GR) e ra 
(68) 
RREA RAE {Simpson》 公 式 ! 现在 我 们 知道 了 用 辛 下 


生 公 式 求 定 积分 的 近似 信 时 ， 误 某 不 低 于 A 阶 ， 即 二 + 阶 。 可 见 
辛 卜 生 公 式 比 梯形 公式 精 信 得 多 1 
我 们 从 梯形 公式 出 发 ， 外 浴 一 次 就 得 到 精确 度 高 得 多 的 辛 P 
^E. Aeub5RHETIIAaU Re E REA. PI, dy (66) AnA 
T=7 AC AC htt OU) 


= ENS z de 4p a 
F OA + T6 OCR 


= 4L 1 z l 4 B 
nacen 15 i^ 十 p58 ci^ 十 OCA?) 
于 是 
af. bi, TCI. 
? 4 
-(brcH sla ROMS MIT 
' 4 ']16/": 
zo 208 * |] g 
! vo edet agg) 00 
4 
b c 
| Ell 00 —— ax 
人 
0 UA 
(9 tiel 3g; - 9 
得 
| TEN 
54d 
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TH 

c 5c A5 .. 

= uox Claras C.f b 

ai^ 16 i 十 ps Ba 了 十 加 [站 
所 以 

1 4 54 864 

LE TEM SOHESS T cs 
45! ， g fà 15 Ps SEL (c B. 


l 4 54 1 
本本 一 ees e kis TIE : 
t 45 * 9 /s Tn da Fa 15 20 a t a 


(69) 
来 求 定 积分 了 的 近 做 值 ， 则 
pt =/ + CR) 
EHAR TA Ere 如 果 把 


1 1 
hetas fof tof. Bin ) 


代入 《69) 上 成， 可 得 出 比 辛 目 生 公式 更 精确 的 访 尔 (Boole) ZA 
式 。 但 是 这 公式 太 繁 天 ， 用 起 来 不 上 方 使 ， 在 这 里 就 不 写 出 来 了 ， 

从 以 上 狗 述 中 可 以 看 虫 ， 将 近似 算法 进行 “外 推 处 理 ? ， 往 
往 得 出 更 精确 的 公式 ， 这 也 算是 一 种 “ 精 加 工 ” 芭 ! 在 数值 计算 
的 评 多 方面 ， 用 外 推 法 常常 使 我 们 找到 更 好 的 算法 。 外 推 法 基 产 
生 新 公式 的 一 个 泽 泉 。 


32 外 推 法 在 常 微 分 方程 数值 解 中 的 应 用 


在 本 节 中 ， 我 们 将 说 明 怎 衬 把 用 差分 法 得 到 的 常 微分 方程 数 
值 解 进行 “外 推 处 理 ”， 从 而 大 大 提高 精 山 度 。 

1. 初 值 问题 的 外 推 

对 于 如 下 常 微分 方程 初 信 问题 


hun DER 


(4) 
yD =y; 
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AUFS: MREFA O, 在 区 域 
D. ugra LUYUTTD 


dukib, Hide 9 有 界 ， 邵 有 型 之 9， 使 


SE uy) | ond 
fEDE 5E D pm3r, Mai (A4 Er ESLa, 5] 上 有 了 唯一 的 连续 可 
Qi SEO, 
WE, BJr, y) 在 区 域 
D, O&K, —ooxg«20 


中 连续 且 偏 导数 3 人 -有 界 。 求 A4) dese [0,12 的 数值 解 。 
CO KAARS (O MEUSE. WO, IRn 
个 相等 长 的 小 区 间 ， 它 们 的 长 二 -…。 分 点 为 


ti 
z=, T= A, gd e. i pase 1 


因为 


所 以 方程 
y' (r= jir, ylr)) (yC0) 2 9m) 
可 用 方程 


Et OA) js uy) 
CA") | 


CH= 5,25, (n — AD AIT 
yC0)— us 


近似 地 代替 。 
为 了 方便 ， 将 点 节 一 0 站 27 5 一 1 的 集合 记 为 
Q,—ir|r-0,5,2h,-,(0—1)A 
但 《A 的 解 应 当 与 原 方 程 的 解 有 区 别 ， 将 它 记 为 yz) 比 较 合 


© £$32315toerBR.BilirschWr &Introduction to Numerical 
Analysis) 3b 7 & 13 
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这。 于 是 得 到 差分 方程 
A rci), 
yg, (pU, 
(25 解 差 分 方程 〈《.4) ， 求 出 wz) CEND, EE TE 
为 问题 CAD 的 近似 解 。 
HE GD 很 容易 解 出 。 记 
YK A) —Y y (C —0,1,2- 0) 


CAD 


Hn 
人 
i C40 "pA 
-了 一 wp 
FieYuqüfto,F.) 
| Y=Y,+ fk YY 
1 dm 2 RfCGA,Y,) 


Y eum, TAFfCK h, Y.) 


y = EES Oen Fal 
这 是 一 种 递 推算 法 ， 一 下 子 就 把 了 , ,了 了,,"… ,了 ,全 解 出 来 了 。 由 
此 看 出 方程 CAO 的 解 存 在 且 唯 一 。 
在 理论 上 ， 可 以 证 明 : 
EIl, VERR D asesb, -oly Lte, AN 
连续 有 界 偏 导数 ， 则 方程 
-y+ R=y (rir yr)) 


: rcÓ,—iz|r-a--iA,12—0,1,--,(10—1)5) 
o | OD (5 -5—) 
的 解 y,(z) 有 如 下 的 浙 近 展开 式 
y, x) rt AW rt RUF. (rt 
TANW QOO) t AY E.G, 5) (70) 


TE JE SX rp, yir Eh fü 


Aai hiphiceducovippbs 
* 工本 。 

人 stg) 

yCaym ys 
ÉIER IIO) HALO) Se LIP CO HEIS ATA BR BRL 数 。 它 们 
HE Le, b) bdüokX HYS. X 

W etag K =1,2, e, N 
RWE, Cr AORE 
Tp 


[Enn GE M Go) (ERIHEH h = h, TL 


T 
n -1,2,- m3) 


(70) 式 给 外 推 靶 提供 了 理论 基部 。 
C35 现在 可 以 进行 外 准 了 。 由 (707》 式 可 知 ， 若 Fo 
ERD 
D, OSTI, —oocy« oo 
中 有 4 NxESUÉg RE N= 20), WEH FE CAO 的 
fiu Cr YE BEIEEJEX 
p, r)— yr) o A x) AW x) A E aA) 


— y( x4 AW (x) - OCA! )8 CT1) 
于 是 

ValT)= yle) + IV (x}+ OCh’) 
Ht 以 

2ya( 2) - uu (x) Hr) OCRT) 
E 


y'Gry-2ysCr) 7 y. GO) (72) 


© $381. Stos ER. BultirschBr & Introduction to Numerical 
ÁAnalvsis PS To a db 
Q vo Fainga biS c. HE. XE. x0 dECa.b OE PUPA 


* d 
| 了 zl ta Eal <H CRO 
AW aim + a* FE rh = ON 
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Dx 
y*(r)e gx) QCA*) 
H GD RAR 
y, vr) yCr)- OLA) 
所 以 ， 由 外 推 得 到 的 近似 解 y*c CAO BIENES IBS TI 
举 一 个 例子 。 对 于 初 值 问题 
Wan TENET. rC€í(2,1) 
y(0)-1 
求 它 的 数值 解 。 
第 -- 步 ， 列 出 差分 方程 
ri yr thy r) (r=0,h.28, (nz 12A) 
| (^ =>) 


fi 
再 解 此 方程 (A,)。 令 
Y à cy RAD KR 0.1, un 


CAJ 


CA | 


将 (A41) 表 为 
人 (Ko=0,i,2, ,nn 1)) 
"S 
于 是 
S 
y.—GcriY,.-—-O-c£) 
Y.—(1ràÀ)Y,-QqQcrA)* 
x 
Yp p= +A)" 
Fy 
Ypa Y aA) 
所 以 
Y g=(1+A)* (K =0,1,;, 有 £73) 
这 就 解 遇 了 方程 (41)。 


现在 进行 数值 计算 。 取 
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n=8, $= 1 
E 
由 公式 (73)， 得 差分 方程 (41) 的 解 了 kx 一 (1 二 过) 
Es. 
Y 1.125 


| Y ,-—1.265525 
| Y ,—1.4238281 
Y ,—1.6018086 
| V .—1.8020825 
| Y ,—2,0272885 
; Y,—2.2806973 
t Y,—2,5657845 


38.41 8. 5m CAD p sé ue RR Rb 
yir) =g" 


A3 - 0, us A A Lb opi: sl. 求 出 
y(x)-—e" TEXAS IS px Ab i6 dep T 
BL 
| 4171.1331485 
| y, 1.2840254 
ys=1.4549914 
ys=1.6487213 
ys=1.868246 
ys 2,1169 
y=2.3988753 
ys =2.7182818 
对 遇 这 些 数值 ， 可 见 差 分 解 yr 与 精确 解 yr 之 差 四 。 
为 了 外 推 ， 再 求 出 ys(T)。 由 《73) AMIR 


© 在 计算 jx 有 时， 还 有 由 区 多 五 作证 成 的 误 痿 。 在 本 书 中 , 藉 不 守 论 这 种 "4A 
R2” Ja 
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yii (1 — 09 


gu mm yalph) €p70,1,2,77,8 Ad gu UET. $ 


SOT ELIT 
ys(2K | 1+16) = Yr 


于 是 


出 此 得 


;4 一 1,1289062 
:一 二 .2744293 
=1.4387111 
4 一 1.624170 


,8335355 


acu 

a m 
Il 
- 


LJ + 
rab meje j= wf 


| 
co 


0698808 


li 
to 


.8367114 


a 
i 
L2 


= 2,637928 


第 二 步 ， 将 上 面 求 得 的 两 组 差分 解 进行 外 推 处 理 。 我 们 用 公 
式 《72)》 HE. dB CO X3 
k —72Y kas Y xà 
在 这 里 
Y kom Uu CE A, V kac Y GER) 


由 此 得 
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|yi22x1-1-1 

Y |::2x 1.1289062 — 1.125 —1.1328124 
Vi -—2»x1.2714203—1.265625-- 1,2832336 
11-2» 1.4387111 —1.4238281— t. 4535941 
i Y*—2x1.624170 — 1.6018066— 1.6465334 
, V 2252x1.8335356 —1.8020325 = 1. 8650387 
| Y $—2x2.0698898 —2.0272855.5 2. 1124931 
y 1—2x2.3387114 —2. 2806973 — 2. 3927255 
Vi—2x2.637928 —2.5857845—2.7100715 


Bé LSRIBBUGSAZRAY aY Qu. ERYS m Hg Hs 
y(Gx)- eB PEE, MEER URNA (Hed A 位 ) 


$8 


| 

0 1.0000 1.0006 — | — 1.9000 | 1.0000 
L/B 1.1250 141285 i 1,1828 | 1.1331 
1/4 1.2556 142744 | 1.2832 — | — 1.2840 
3/8 1.4228 1.4387 — i 1,4538 E4550 
FE: 148018 1.6242. — ! 1.f465 1.5487 
af8 1.8020 1.8335 i | ,BB850 1,8633 
3/1 2,0273 2.0889 | — 23.1125 2.1168 
7/8 2.2807 2.3387 | 2,8927 2,3983 
2.5658 2.68375 | 2.7101 2.7183 


dENEXM E, 4$e41X GÓC SD/MEHRUIS 大 优越 性 。 我 们 把 误差 较 
大 的 差分 解 了 *, 和 了 kew 用 外 推 公式 〈74) 处 理 一 下 ， 只 花费 了 
很 少 的 劳动 ， 却 得 到 了 精确 得 多 的 结果 1 

2. 进 值 问题 的 外 推 

对 于 如 下 的 常 微分 方程 边 信和 问题 


—y'(r)-ag(r)—f(x) rC(0,1) 
yO =y yD =y aX Raz o)D 


D 4e> 0, jin Æw, HEER D 的 解 存 在 ， 唯 一 。 


| 
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(这 是 线 狂 常 微分 方程 的 两 点 边 值 问题 中 
我 们 打算 用 外 扒 法 求 它 的 数 什 解 。 仍 到 三 个 步 又 
C15 FRS CB) 相应 的 差分 方程 。 


将 区 间 0，1] "等 分 ， 令 请 二 ， 分 点 为 


| 
HH BIETER 
CyCr,,,)—2pC€rz 0 o nGr; 0/8? 
AXE AUAM LE E Ay (r0, GAANA 
人 
(GI) (一 一 age) 二 2 
|: yu COO — ya, Yay, CI —1,2,-",n—1)5 

(25 fu JJ ER (BO ,oRdUEm Su Cru), B E TEN 
问题 CB» 的 近似 解 。 

Ji fid (B) 是 比 性 的 ， 它 的 种 存在 且 唯 一 ， 这 一 点 只 要 
观察 它 的 系数 矩阵 就 可 以 肯定 ， 不 详细 说 了 。 

在 理论 上 上 可 以 证 明 ， 若 【已 ) 力 程 右 端 男 数 天 >) 在 [0，173 
X8 4BpDESESEEEID, B CAD BUSCO ETE IBESEXUSOT 
X 

gx) yCG) VW xt h5Eh) (75) 
EEA, MIOr)ERiS ACE, CEDE 4 Epi ROG, PEE 
五 :rzr) 有 界 。 我 们 在 本 节 的 最 后 来 证 明 C7) 式 。 初 学 者 和 只 X 
心 外 推 法 的 应 用 的 人 ， 梧 以 不 看 这 个 迹 明 。 

(35) 5E. ET OD 这 个 斯 和 近 展 开 式 ， 就 可 以 外 推 了 。 
和 过 去 的 做 法 一 样 。 因为 

h 


hy PEE 
yalr)= yte) + (5) Fie + (5) E( 2, 2) 
O 方 各 由 y" (+f) 的 系数 取 为 ~1， 是 现代 数学 工作 者 为 讨 沦 一 些 再 论 问题 Cun 


< 广电 解 ? 问题》 的 方便 而 设 的 。 和 如 果 您 不 导 慨 这 一 护 ， 可 将 此 方程 的 每 一 
上 项 光臣 上 以- 工 。 
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—y(r)d FW) 十 OCA) 


X C5 式 可 写成 

gx) yr) RTI x) 4 OCR*) 
若 取 

g*(z)m y CAys Gr) — G2] (76) 
pli 

y*(x)-g(xr)4- OCAS) 
还 是 举 一 个 例子 。 对 子 边 值 问题 

—y"'(zx)ty(r)—-0 r€(0,1) 
y(0)20, y(1)=1 


用 外 推 靶 求 数值 解 。 
第 一 步 ， 写 出 与 《8》 相 应 的 差分 方程 。 


将 Co, 12 区 间 # 等 分 ， 令 有 二 小， 分 点 是 


n* 


a» f 


Xg— ü yt = fi Q24,— Rh. ax mnh-] 


差分 方程 取 为 
(B) no d o 0 dor clic c 
y,(0)—0, y(t1)=1 (j 21,2,-:,n—1) 
也 就 是 


(B,) pa Mola A ua Gr) EYLE) 0 
y, C0) — 0, g, C121 (Jj =1,2, "=, —- 1) 
第 二 步 ， 解 方程 组 (3,) 


取 # 二 2， 则 有 二 亡 。 这 时 只 有 三 个 分 点 ， 即 


c 
YCE = 0, y,(x;)—py,(10-—1 
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RA MERETE EEST 
E (BD h, Risi he. 


UMCT m 2 £l)yGossGo- 0 
Bru 
- Tyr)+1=0 


Wn m 
于 是 


^h TUE, ER vr) RK, M 


n" 4, h’ =-7 


将 区 间 CO, 12 4 等 分 ， 分 点 有 5 个 ， 即 

Xu 二 0， zi 一 地， r=, n-l, T= 1 
ia 
A 


2 
MAEA CB. EA 
Wp 
Hoc 0, y.—1 
详细 地 写 出 ， 就 是 


1 
c EL 
h= 1 


E y,r(ir;)-—Hi 
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于 是 


y 
| y.-0.6999662 = w.( 
第 三 步 ， 外 推 。 用 CO ct, fbr. 处 外 推 ， 得 


oh) 


1 4 
ME (4 x 0.443680402 -过 | = 0, 4434257 


SEXIES F. BIS CAO 的 精确 解 是 
yle)=zr lete) 
在 = 过 处 ，y(x) 的 值 是 
sat ei— eT) 0.443409442 
外 推 的 效果 是 明显 的 。 
由 76》 式 可 知 ， 用 这 种 方法 只 能 在 点 x 一 十 处 作 外 推 《在 


其 余 点 处 因数 据 不 足 ， 不 好 外 推 ) 。 这 就 使 数值 解 在 点 z 一 十 及 


= = 守 处 误差 较 大 ， 造 成 数值 解 在 各 点 处 攻 差 不 同 阶 。 最 近 , 我 国 


n 


[ WI TM 
" 
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PARERA EELT- RAMINA, MATAN H 


UL 3k dt ABENE. HERRN ARAARA LAE 
在 分 点 < 一， 了 处 进行 外 推 ， 
用 分 点 


.ns 一心， ys ~z, =nė =] 


PRE Co, 1J 5 等 分 BHIA 


h n h 
£ g = Ü, aT Toe" os Tei 
h 
= ne, = i 


E [0, 12 2n'$4r, din ^52 BJ EISE i 
£2. PIU ML 


把 22 等 分 的 分 点 的 混合 记 头 
a=] ME Nn !?- 0,1. , 2n} 
如果 


Ws ZT， 
但 是 zi 和 TY 的 中 点 
T, 二 之， *raQ44)/2€ aN 2D 


XxxHJ, 32114.254) 7: REC B PR RE LH 
DIG iba) 
和 
ysCr;), BC. Yal Erg) 
这 几 个 值 。 于 是 ， 用 (76》 A+hiE, [Du 
ysGr) mu Gn) ) 


j 
yir) =} 
m» 


Ò HtA AAT tief RE. JEZGRCEZpmipirsk 893m. 
Q MAND AEREA Qin SETO ELE. BAr 属于 Qr ERAT 
Bus 
Gn Heck TeMi aAa 
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IN NTNES oo PVID. Nt edu.cNnAVIDDDS, 


yen d= 4y&Gn a)r) | 


但 是 ”g(x,) 的 值 不 知道 ， Br BEBE JH CIO) AEQ hsh HE, DA 


£6 385 BOE AS s 
ytisom yii e Go TG EAER, 
6 " [3 
ml ct ) ] (77) 


现在 用 《77) 式 在 分 点 xz 一半 及 + = SEINE, d 


si) vi) T 5i y 4&0) — y, COD 


—0.21511777 + E (0. 443680402 - 3] 
—0.2149903 
TES = iit (1 — (iy i 
(i) TE 3 C3. Vil n» 
7.) 0.699662 ü st 0.443680402 — 3) 


= 0,6998388 


在 此 二 点 处 ， d dE u(r)———— (er 一 es) 的 信和 是 


e*— 1 


1 e a E 
—.d————s lo i — ý mE 
(1) zí E€ e €)-:0.2149523 


3 e 2 Zo 997 
-= 一 i 一 4 z 
(1) E rte e 43—0.6997241 


将 以 上 数值 列表 以 便 对 比 
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444 
2 
z, i E EN | 
4 i 2 d 
TREE | | —— | 
gs Cru) | 0.2151177;7 | g.44365DAD2 | r,5989552 
2 
b rx 20.21 45303 0.434134225T7 6,C955388 
YE) 0.,21435723 0.,4434024472 O.tUuS9T24I 


M. X& 4r 91 3-783 DUE E E TES 08 23 23 ER P DA E GE. 
HER COD 式 结 出 证 明 ， 
对 于 问题 


(B) se r€(0,1) 

YCD) = Ye YI) =Y: a 为 常数 。 现 在 设 2 之 0. 
这 是 线性 常 微 分 方程 进 值 问 题 。 在 微分 方程 理论 中 己 经 证 明 ， 右 
fG)fELO, 13 连续 , MCE) 在 [50,1J 中 有 了 唯一 解 y(x+)}。 这 时 因 
为 i 
y'(x)—ag(x)—fCr) 


mylo nE, Mu ett, HETE, #7) 
yP lr)=ay" lr) f" ix) 
所 以 zz) 连续 。 又 
O(n (2) 
H Wyt zz) 也 连续 。 这 就 是 说 ， FO OPTED CT) 
MEZI rela, dD., Rip BOX ME-—H EVO) CE 
在 [a , b M8 6 Brie Sk 5j 3C. 
现在 ， 证 /zz) 确 实 具有 4 阶 和 连续 导数 《 当 zErLa,p]) ， 我 
(IX urHj (75) 3X, BP 
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yu Gr) yla yA ATTE Co 4 EC. Rh) 
XI CIL EY) o IE a) 

OHET-PB Tj, JPA OMIJE r) 
EEP, vO ECS OE CI ot , yere) Bj f 
CEDAS TAR s E CAT 9t, 

EL SIFE PRH o Sg». HD a,b. 
l4 EDEESE TE. GU TSOGAR (XV AQO HU BRAJ fE 

eiu Pu gj—wmx: eX e. ax 


tay, us de 


id 
二 
TA 
—Py(r.)—h* IF: 3 ,) P Ea) duda er 
(782 


füyCryf 6 EE H 


Å 


ki 3 
UE) = hy {rr T yir, +t a ue (r;) 
E 


SNR AË ns T 
ty MiB jt--y"ceydpB os 
41 5! 
| " " ^p At, di 3) 
YLT) = yt) —hy (x, "EA (x) — ay (m,) 


TT A? 8) 8 
ty rt OC) 
4t üt 
于 是 
Pyr = Tyl; 29x yle, Nk 


= y'(r;)-c hy (r) +O) 
IWORA a 阶 连 续 导 数 ， 所 以 
FG ja) m Gr) B^ Qr) e ER C) 


a 
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127 。 
a AE 38) ( p 1) -- 0C A) 
W (x; mW ERETI Cer) + oh Gao) 
-g n" 0 (x) + OCHE) 


-已 A. 
PA 


PLx SCH (n, 7 21V Cx ,) -W (m, 91/7 
— ga AEG) t OCA*) 


—H'(xroco(h?) 
所 忆 C78) 化 为 


E O2 ry CAS — B* CHI" Cr OCR YO 


127 
—A PEG, tayia, hit (c, ) c ah* 
x Bg reris. 


B 
— y" lr, )— fiz, )+aylz:)+ di - bu o (x) 


NC LUE IE hi0) - P E,G,) 


+ab, izr, 0s 0D C783" 


ihiksp us, GWIVECO mE 


1 

一 于 7 | nr ere rd) 
(C) | Cro atl Qa 12 Cr) 
‘(0)= 0 


Jie) ABpESESEX 又 
(LY -PEE (Cfo=1l1.2 499—345 
E U= 0 S= El) 


T aaco ta 

2, TYT yi) +H igi + 
] Coa N= By amg) SH n = 0 
Dn 


SE. Q0 H. Siu, Salh 


Jy. im =y, y,tn-gy'b 


oat r+ AB, E G an 
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则 (78)' 式 成 立 。 在 这 里 ， 注 意 
—y'(x)crtanyg(z)-—fíz) 
Br 以 
—-g'iz;)tayg(Cz)ef(z,) 
自然 成 立 。 
观察 方程 (CC》， 这 是 洋 桑 数 线 性 第 微分 方程 。 方程 右 端 的 
CLA A Kcu (^ COR 2 MERER DUERARI GOHA, E 
由 方程 CC》 可知 WV(z) 有 4 阶 连 续 导 线 。 
再 观察 方 程 (DY ， 羞 它 是 否认 大兴 和解 让 ,Cz7)。 把 ( 品 ) 详 细 
ep Ex AE 2) b Fr A def renti 


Bil 
EAr 2ERx Ea a E uCe;)y-00)4t 

—0CA?) 

令 

7 Bp pem 
网 
=E | 十 有 和 一下， 十 下 下 天 ,一 站 下 2 C5 

现在 考察 如 下 疗程 
ao ud" ale er DN CED 
E =Ë, 671 — 1.2 ,,8—1) 

B 
(LaSi E t ou gau ce. 

X 


"nl 4 — 


AUGE oae E E, - Ean E, E, t Ee., 
J=] r=] 


a~l 
zd c UB pe E p Pen E ELIT 


了 到 二 )-1 
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* 129- 
r-? aal 
=Y (E; EnEn -5 (E, -En dDEr 
i= ù nel 
= X (E; En Em Ent (Eos E)E 
j=l 
—-(E,,—EJE,A 
=- SE Em) +E EE- Eo) 
—(GE, ,— EC(GE,  ,—- E 
-— VE, EQO GERE.—0 一 也 
jn 
Br 以 
并 一 了 二 一 F- } 
SEQ Em) + Doke DFE, 
Pe j—1 i=l 
于 是 
XE ES +Y ak E= DPE 
je i=j 
X 
j-! 
|E; =| BSE, — E| 
D 


n= i [PM 
eX En- EL «vx T SG. EO 


(D XHJmERSCauchy E ab. xY ya a: Y ybt 
i 二 1 


i=l 
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Erb 
a~l 4—] n= 
HE j| tn S OE RS RES n- ak pN - Sr,E,) 
1 一 由 1 =0 i=ù 
于 是 
i=l tUm. 75 E] 
£o EJ S uo M M JUR 
lE | Sn OE > er 
«n V nmax ip V n max |E, 
Í i 
(H 7x n= 1) 
EB! 
max IE] *2en'max q. "max [LE 
所 以 
(max |E )'sm'max (9,| «max | 如 
i i i 
由 此 得 
E 
pcd El s E max Ig; C793) 
5:8 n- ， 


《注意  nh—1) 
(79) 式 是 关于 方程 组 〈 五 )〗 mHEEBRI— IRA. GRECE 
方程 组 (HD 的 解 存 在 且 唯 一 ， 只 要 观察 《〈 瑟 ) 的 系数 矩阵 就 可 
在 方程 (E) mp, RP, =A, gb RJPECOD'), i O9 
x, n ACD E—ERE;:-—1,.2.7,5—1)872, 


振动 论坛 http://vib.hit.edu.cn/vibbbs 


Mu, M0 OMA. o) hdd OSAROH k 
20) Hj, NE 


max jE «M 
05 rac -] 
HHE, = EIOS. 于 是 LTO) 式 得 证 。 


$3. 外 推 法 在 仿 微分 方程 数值 解 中 的 应 用 


在 上 一 节 中 介绍 了 和 常 微分 方程 数值 解 的 外 推算 法 ， 这 方面 的 
工作 已 经 相当 成 处 ， 有 了 许多 好 结果 ， 在 国际 上 也 早 就 为 计算 数 
学 家 熟知 了 。 人 了 但 是 ， 如 和 何 用 外 推 法 来 求 侦 微 分 方程 的 数值 解 呢 ? 
在 国内 的 载 科 书 中 还 没有 谈 到 过 这 个 问题 。 国 际 上 已 有 一 批 成 
果 ， 目 前 中 、 外 数学 家 还 在 紧张 地 研究 者 侦 微 分 方程 中 的 外 推 
法 。 这 个 领域 中 的 外 推算 法 ， 在 理论 上 难度 较 大 ， 在 实际 计算 方 
茄 也 比较 复杂 。 在 本 节 中 ， 只 对 一 个 具体 的 偏 微分 方程 进行 外 推 
计算 ， 并 且 结 合 这 个 具体 问题 介绍 一 点 现代 绪 果 ， 使 读 着 初步 了 
解 偏 微分 方程 外 推算 法 的 一 些 特点 。 

对 于 


一 二 了 TV (XV)E QUO 3E TRI E RIJESE 
m R JE TR 
uj aa — 


Ei, Au = unt am- 是 二 维 Laplace 算 符 )， 符 号 00 表示 区 
RIHAR. fOr,y) J& TE DOE Q ptr gg ^ 169 CLE PE, 问题 14) 
就 是 ， 求 孙 数 u(r ,y)， 净 它 共有 二 阶 连续 偏 导 数 ， 且 在 盾 形 区 
域 虽 内 浦 足 方程 一 -Ju 二 fz,y)， 在 矩形 的 边界 9 包 上 满足 边 值 
Pelr, yas 0 CRER l EMERE = QUOOE 
连续 ) 。 方 程 


CAD 


—zZiu(r,y)-f(x,y) 
是 天 家 熟知 的 Poisson 方 程 ,问题 (4) 称 为 Poisson 方 程 的 第 一 边 
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Ei RURAR OO R JHT 
Je HH SPIESLTES EB je Ria 
现 人 在， 将 问题 《4) 上 其 体 化 为 


—ziu-—zsircosy WMir,y)c£40 


CA) l | 
ul 15 0 Q, Ox 《了 在 名 是 正方 
JE " 
CSS, E. 
第 一 步 ， 写 出 与 ( 4) 档 应 的 差分 方程 ， 并 求 出 此 差分 方程 的 


解 。 
HE o ,区 间 分 成 个 相等 的 小 区 间 ， 分 点 证 


2 二 


又 将 [ 一 六， Emme 4 等 分 ， 分 点 是 


JE Tọ JU 
Vu 3? Qi——t hs Wie TAB quus 


-— A T. a- 
A t 2 ( 3 
SIE IE; E QA WIESE S. CHILA R "d 
A? m8" A43] 
(x,,9,) 15,525 0,1,2 "-,n 《其 有 (二 1 个 节点 > 
id 
ur ju 
(L7 (a SP RC eU XX 
ur, y, oe H jana uus, H xr; a.a) 
esu, — uu) 
X, Cp da a A E TUL XX 
Us WH EI Ya) uK Ey NA 
AENEAN m Cia m isa fA 


z 
=| Hpi T 2444 T PE 


HAE, A 
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Uyy EEY SCH th 
所 以 
Zfuj,— E yu) usum gu) 
= reru Micquc Wut Mac du aA 
(80) 


E (CO) xo XR ARE OG Yy ORE E B. AX. Gg 
T AALA KRE XX cA. COPA AARAA” o du 
Hiri t Hant "asa E Epeei TENi OU, 
jgBod3—:3-SHF. 5 CAD PERAI 2:542) Jr fan] ELS x 
(uu Q,kR—1,2,-—,-1)0 
u;,— 0 CO k= AA e 0 3X n 


现在 来 解 方程 ( 4,)。 先 到” = 2， nh, 这 时 


CAD 


Om 
t= OQ, TQ——, Ff 


2 
= 一 0,90 
Ua 2" 1 ”Ys 


X, 
OH S Hp Wort Win Hs CO— AN, 


FE D 化 为 
feit at uin t TEE » sin 7; 0080 
hil eed) m eH da 

实际 上 上 只有:# 是 未 知 旺 。 册 此 方程 解 册 xi， 得 


a 
TN 一 车 二 1.2337 


方程 CAD 的 精 衫 解 是 
u(r,y)-—smnzcosp 《这 可 以 直接 验证 ) 
出 此 得 


T = 4T 
4C. n ) — Sn 一 co80 一 1 
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Tt T VEN 
"t RA (Z 0 )o soci 22278 (7 ` 0 ) = 1 的 近似 值 ， 
误差 太 大 了 。 怎 么 办 呢 ? | 
我 们 很 容易 想到 一 个 办 法 ， 就 是 在 求 出 #( 态 ，0 ) 之 后 ， 尘 


求 出 wu. o). mmm 


sb m o) o) et o) 


进行 外 推 。 也 许 这 样 可 以 得 到 更 好 的 近 伺 值 。 
虽然 我 们 现在 还 没有 把 担 ， 但 不 妨 一 试 。 反 正 间 题 ( ABUS 


确 解 已 经 知道 ， 可 以 直接 看 出 这 样 算出 的 4*( 瑟 ，0 ) 的 误差 ， 比 


si 号 ，0 ) 及 (Go o ) 的 误差 究 况 是 大 些 还 是 小 些 ? 这 是 一 
个 有 趣 的 试验 。 数 学 定理 往往 是 通过 大 量 观 察 才 发 现 的 。 现 在 ， 
我 们 来 求 xi 号，0 )。 为 此 ， 写 出 步 长 为 "= 方 的 差分 方程 
VENUES 


HR HR inti rt A O22) 
Ar aa. | 

ví cms. Quien) 
这 相当 于 将 正方 形 名 等 分 为 本 二 16 个 小 正方 形 ( 取 = 40, A 


以 按 下 式 夺 代 求 解 。 
| 
| 十 at E "- y ) T ( "est KJ] 
J a sin (3 了 ih 
upo OKZ- AoA Q,K-—1, 
a 3 3 ) 


AdüufDoupGmOg.E—1.2.3) 
H jef ko u TR i SGK 1.2.3 BELEIRE, LI 


振动 论坛 http://vib.hit.edu.cn/vibbbs 
-135+ 


u Pe EAG: BEE., D 
现在 ， 取 :一 0.0001， 得 到 ( 41 的 解 如 下 ， 
Hi =0.5265 5s,,—0.7446 wi,s=0.5265 
aa 一 0.7446 ay 一 1.0530 ty 一 0.7446 
ra 一 0.5265 dss=0.7446 as 一 0.5265 
于 是 
Ze 0-5... 1.0580 


现在 采取 第 二 个 步骤， 试用 公式 
iC) yen o) 
CL. o ERNE, 8 


(C, 0) C4 x 1.0530 —1.2837]—0.9928 
SEES (Z o) = sin 于 cos0= 1 比较, XI (I. o) 


uL. oy “区 王 ， 0 ) 的 误差 都 小 多 了 | 
这 种 算法 有 什么 理论 根据 呢 ? 许多 数学 家 探讨 了 这 个 问题 。 
苏联 学 举证 明了 如 下 重要 结论 ， 
MEHE D 中 方程 右 端的 函数 Fr ,ob 满足 
Fir, yE c?**(Q) Coa 1) 
H.EQBS UU fA, A 
fix,y)— off. .—f,—0 
则 问题 (<4) 有 了 叭 一 解 x(z, y)€ c7(Q),. BL 3; RE EHCA B9 IR 


uC, yA ATRI 
ulr, y= uir, y) +t A'W ir, y}+t E,(r,y) 


TRAH 


T YRAETA AAEE Aak. MiA, SàckRADEDTOE 
A Hp PEA EIE. Aie it RALEA YA Ga uss Ai nt US RA 
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IE. gy sicht, Wir, yE cY) 
这 是 用 简练 的 现代 数学 语言 表述 的 结果 。 在 这 里 作 一 解释 。 
介 是 条 形 开 区 域 介 加 上 它 的 边界 3。 所 谓 
fiz y) Ect Q) 
EHS vC B — ET DRAR RT RE LERT, 
HF RHEA SE a o r y ET yf as faer fays 
fa PRHE — AO EJ 
PCr y —9Ga uii EME (G7 z, Hu Y 
(zi ,#1) 及 (xX,,ys) 是 人 中 的 任意 两 点 ， 必 是 某 个 正 数 ， 
u(x, y) Ect) 
的 会 义 可 以 类 推 。 
这 些 条 性 看 起 来 比较 难 掌 握 。 但 是 ， 我 们 指出 ， 上 只 要 
f(r, Ec NR) 
GEG, ERLEEN A 阶 的 连续 偏 导数 ) 则 
flr yp Ec 人] 
所 以 ， 我 们 可 以 写 出 一 个 出 较 实 用 的 结论 ， 即 
Zee fai, 
(10. fGr,y)€ c*CO) 
(2) TER P 0 Bh V3 f e b, fix.y)- 0 
(32) EEEO RA, fef 0 


al 
C1) BHBECAYdWE—fu(r, y)€ c9**(Q) 
C225 HÆ Rs, DARRERA 
u,(r,y)—u(r,yg)ct R(T y) tO) (81) 
EKER HHF r, yE 


W(x,y)€c*(Q) 
XARES HENE. SEE B. d GD x, $$ 
Q "fci EMO LUSEUGX—BUEITUNOEXGYEXESCRGRRT. Utd WWE 


HA’ dm)GilbagE Trudinger N Y Eri Sd RB y 第 一 
Ger did e nk RAM E € ga gR a on Eig, 
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hM T 
utr =H y+ |5 Hr ug) 4 0CA5*7) 


Pr VJ 
d 、 1 EN x x 34-4 
u*(a, yi lelas g)—H,Cr. SES ulz, yt 0CAT*7) 


Hik F., RECS FEA OPACARA yA E E 
BRG (DRG) MIT CAL 00 22 2 Eu Mua Cm. y 


MRAR ua, y) =A Cu Gn) cu G 2 (82) 
2 


Hribi, miu", yW E 
lu*(z.p) ule, y) 20Ch?*?) 


在 刚才 的 例子 中 ，z ,一 2sinzcosy， 它 有 任意 阶 连 续 偏 


J 


S3, Xf. fafaa, o. -—). (o. z) (s.-2) 


a(r, TAP, WRA ORERE), COAC), 
所 以 可 用 《82) 式 作 外 推 。 
MAREA o ) 处 进行 外 推 。 在 其 余 点 处 因 缺 少 x,(z， 


DRE MUREA (82》 式 作 外 推 。 但是， 在 这 些 点 上 ， 贡 可 
MERARI. CO 式 是 关于 一 元 通 数 的 陈 - 穆 公式， 
对 于 二 元 了 盘 数 有 同 祥 公式 ， 在 此 就 不 详细 说 了 。 
REZE- E OFAR, y WEE EEA IE 
(1), (2903), EEP (820 式 进 行 外 推 ， 可 能 越 推 误差 越 
Xi 所 以 我 们 在 解决 实际 问题 时 ， 决 不 能 不 管 条 件 任 意外 推 。 
3419 XR T XUEXE ERI Poissongr ERE TEE. 实际 上 讲 
的 还 是 这 类 问题 中 的 一 个 特例 ， 即 正方 形 域 上 的 零 边 值 问题 。 对 
于 更 复杂 的 区 域 ， 它 的 外 推 靶 也 讽 复 区 ， 而 且 还 在 不 断 改 进 。 关 
于 “ 偏 微 外 推 ? 这 个 故地 ， 现 在 好 像 还 挂车 “施工 重地 、 亲 人 免 
进 的 牌子 ， 苏 联 数学 家 马尔 丘 克 和 谢 多 洛 夫 所 车 的 笃 关 分 及 其 外 
推 训 一 书 的 英文 本 383 年 福西 方 出 版 ， 在 此 书 中 对 “ 偏 微 外 推 ” 给 
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d y —ÉÍRÉfBUSNIB. EBAR. Bh At A go b e a p ue] 
未 作 说 明 。 总 之 ，“ 偏 微 外 推 ” 工 作 还 远 远 不 如 * 常 微 外 推 ” 工 
作 做 得 完善 。 希 望 有 兴趣 的 迁 者 能 在 这 个 什 域 中 进行 工作 。 


$4 EJI Romberg) R} 


讲 外 推 法 而 不 提 龙 员 格 算 活 是 不 合适 的 。 

Romberg 在 1955 年 提出 一 个 很 好 的 外 推 秆 法。 用 这 个 算法 可 
以 逐 深 消去 从 和 优 阶 到 高 阶 的 各 误差 项 ， 写 出 一 个 由 近似 值 组 成 的 
* 外 推 表 ”。 在 这 个 算法 中 , 先 算 出 一 些 初步 的 近似 值 ,然后 求 出 
第 一 次 外 推 值 。 将 这 些 人 再 外 扒 一 次 ， 得 到 第 二 次 外 推 值 。 再 将 
这 些 值 外 推 一 次 ， 得 到 第 三 次 外 推 值 。 如 此 继续 下 去 ……。 一 般 
说 来 ， 如 果 先 求 出 # 个 初 近 似 信 ， 可 以 得 到 sn 一 1 个 第 一 次 外 推 
值 ， 再 得 到 ”一 2 个 第 二 次 外 推 值 ，…， 最 后 得 到 一 个 * 一 1 次 外 推 
值 。 最 后 得 到 的 这 个 近 仆 值 ， 它 的 误差 会 相当 小 了 。 

为 了 说 明 这 种 算法 ， 还 是 回 到 求 x 的 近似 值 这 个 问题 上 来 。 

我 们 打算 作 4 次 外 锥 。 为 此 ， 先 要 求 出 5 个 初 近似 什 , 还 是 用 


.odL 
d, = = 
"n 


yalin =d, 6, 8, 12, 24, 18 


. H 
dq? —4sm-- —2.8284271 


， 
dj" —68in —— 3 


， 4 
dq —8SIn'g —3.0614674 
ul 
12 
D 


24 


d?! = 12sin 75 —3.105826 
Ld f^? 一 24sin 


第 一 次 外 推 值 是 


=3, 13262856 
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e () 
(00 4 6 
"T-T OD 
4 6 4 6 


el G) xs 


| gr E'EXE ESTEE 3.061464 "Gy u 3 
| (+) x 5 ) 6 B 
| —3.1404968 
Y 
( 1 ) (Ly 
dv e cl uem .x8.1058286 一 — M27. 


1V Dx 

| e RECO 
1 x3,064674--3.1413171 
i) (5) 
d = AS — x 3,1326285 — — , 

GSH GE 
| 12 24 12 24 
x 3,1058286 —3. 1415618 
第 二 次 外 推荐 是 


| di? ss. G) | x 3.1404968 — ~. (x) . 
ZEN 2 


x 3.1372512—3.1415786 


| e] 
A 6 x2.1413171 = 


= AG Y- 


x 3. 1494968 =3.1415905 


——— € €— 


dr e) , X 3.1415618 一 2s DN 
(3) zt (z) RU 


x 3.1413171— 3. 1415923 
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S8 Ux PME S E - | 
gU 0. —Xx23.1415905 一 (9 


G-G -(à) ire io (2 


x3.1415786 —3.1415918 


"I ICI 


— .1415923 — ; 
3 GYy-() x 3.1415923 G-E 
6 24 6 24 
x 3.1415905 —3.1415923 
第 四 次 外 推 值 是 


Ge G) 


Us me EUMDEM c E 
(人 -二 (1) - G2) 
x 8.1415919 :—3.1415922 
最 后 这 个 近似 傅 qx") ， 已 经 相当 精确 了 山 。 由 此 得 到 关于 的 近 
但 值 的 外 推 表 如 下 


3048475 | 


S374 413171 1 


如 果 您 认为 ar” 还 不 够 精确 ， 那 就 多 算 一 些 初 近 似 值 ， 再 


© më RHH, RS NO i1mBEIESiLXEA E. SS UE Wu 
$800323E, "mM EclRS.1415018, rphdbsrh 2 BN SMEH,. 


* 
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多 外 捧 几 次 。 几 引子 计算 机 进行 外 钼 ， 是 很 方便 的 。 这 种 既 序 很 
eM. xt, SIEGEN Gn Cd. RAE AER p a AGAS 
Xingu— Mm. [EREIIBOEMERNIEXGSIPPRIE BB C 
dxf E TA HP Me JSESDTAEWGSYET5. 

ME EHER byhsRHETRIERDUTÜU;S JE Mi 

A 一 一 h, Wd, sin =d), FË 


dh) — ngin l sins 
7 " 


] ya 1 LAB 
rm T PE [rh P -CHAJ —— »-- 
f 3 HT 


* "r9 
=- ~ ' M ri 
ET NI 
By 
d(Áh- m cft c, hc, B9 eu ht H (825 


XtdRGCOAMEHECRIJPGC. "Ih DS eg A EEROR CAE 的 项 。 
(825 RERI EH VMeDIPUERIB. 
d n 5y SEE EE Haa. tiaa tt, HQ, ttt, (n nue ny) 


TU Ric A, —— B Re SA 


ho hs bu m Vm 


d(h. y, RN at hay e dC urn, 
mikgWikfe—cxBRE. "IO 
d(h Y — mach c,hi-c e, E652 ch te 
di jeter doplA . dug at eaat 
所 E 
h* dh yos mh uon e b RS. roA htt e SAT Au 
Tedtht, Re 
h*d(h,, mh! c, hi h?, + ch? ht, + c, RT RS. 
Tog AT.. ee 
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ACE 
于 是 
AidGCh, ,) — h* dh nh ht eC RRt, 
Rt ca Rif — hh H eA RS, 一 下 8 下 2 


sboss 
由 此 得 到 
hidCR, eohi ou chi ktu ehh? CR? 


i heL, 
BAT. —cQhT AT. Chtt hi.) ee (84) 
现在 令 
ath ee 
[A*dCh, ,,) — At Ah hp, d£? 
Bp 
dí? z[Rh?dq'*P fd ICH? Ag) (85) 
由 (83) 式 可 见 
dq sate h*c ted i eh ee 


所 以 
dé O= a+ h?) (86) 
由 GO 及 (85) 二 式 可 网 
dr =m hkh?) (87) 


这 就 是 培 ，d 7 kedi O REZANIE Z dO HORE SW e ah E 
fü. 

RERA RAA ahid =0,1 ,2), 使 它们 
Hd f"? 的 近似 程度 更 高 。 


把 (80 起 写成 
df'U-—s-chibi. rh hE Ri A.) 
—c,h*h VOR Suum (RR) 
二 是 
dis exc UP. eah? UAE Ad. qp. 
— c.l VB. CRT LIE AT.) —e (80^ 


AEREA KIRRI, MA R (89) 的 各 项 ， HAt, xe 
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(88) HIFA, HEHAD, 15 
Afd pote. de epe tum 
r Vui. CART. — R1) 
sup E S edt opi xq ui cde cc 
c bI oe 5 o o LL 
所 以 


[站 —h3..)—mt— Cah 3 B? Ata 


uus 
令 
gi ad hd WINS AT.) 
fiis uit mo ESI RTI re (90) 
于 是 
di nv ATAT AT.) (91) 
4i BZA OC abii, et (91》 和 GD PA, HD Andi H 
JE H ede 
SPA LUAMEmR. H (O0 A, 可知 
dtt Sar Cy VAT Et corde 
"ur pL 
hid Et mhiad e (GE T Riu) em 
全 


di" ephlg tpi qp A enia) 
dí" =at Olh Al. hllAT..) (92) 
dd 是 第 三 次 外 推 值 ， 
一 般 说 来 ， 令 
dA =Chid i Bnd DA — Aie) {93) 
可 以 证 明 
dEn =at OLA A hi rk) (94) 
(03) 式 告 诉 我 位 如何 从 第 天 -1HE E013 BKK 
外 推 值 d EU. 
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Jü c U5MEfuWpsgede. W ME 

KPRTH S RAME M Milie itii IDfeNPST. 这 种 多 
KAEH MJER r. AATE oi A ANSA M 

HEEE, Bii TEAM EEG, GLAtAAGKcfgrA LEA 
RAR, PAPRATI AARE, 但 是 外 指法 
的 活动 范围 决 不 止 此 。 它 在 积分 方程 、 代 数 方 程 ……，。， fA, 
有 广泛 的 运 岂 。 我 们 在 第 岂 章 中 讲 到 级 右 限 元 法 ， 也 是 一 种 近似 
解法 。 所 以 由 大 限 元 法 得 到 的 解 ， 也 可 以 外 推 。 自 好 和 华 以 米 ， 中 
HAEA RGA NEDÁ, KA TEKI AD, ADDED 
所 限 ， 不 能 谈 这 些 工 作 了 。 

ATAPI EFi APERE HR kR A T 
的 要 求 ， 反 而 刺激 了 计算 方法 的 发 展 。 广 泛 运 用 和 发 展 外 推 靶 ， 
是 当前 的 一 项 迫切 的 任务 。 
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